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SECTION  1 
INTRODUCTION 


Satellite  communications  systems  that  utilize  transionospheric 
propagation  links  may  be  subject  to  severe  performance  degradation  when 
the  ionosphere  is  highly  disturbed  by  high  altitude  nuclear  explosions 
(Arendt  and  Soicher  1964;  King  and  Fleming  1980)  or  by  chemical  releases 
(Davis  et  al.  1974;  Wolcott  et  al .  1978).  During  these  events,  the 
increased  electron  concentrations  and  the  irregular  structure  of  the 
ionization  can  lead  to  intense  Rayleigh  signal  scintillation  at  the  RF 
carrier  frequencies  used  for  communication  links. 

Under  severe  scintillation  conditions,  the  signal  incident  at 
the  receiver  can  vary  randomly  in  amplitude,  phase,  time-of-arri val ,  and 
angle-of-arri val .  If  all  frequency  components  of  the  signal  vary  essen¬ 
tially  identically  with  time,  the  propagation  channel  is  referred  to  as 
nonselective  or  flat  fading.  *:'en  the  scintillations  exhibit  statistical 
decorrelat ion  at  different  frequencies  within  the  signal  bandwidth,  the 
channel  is  referred  to  as  frequency  selective.  Frequency  selective 
scintillations  are  therefore  encountered  when  the  communication  link  band¬ 
width  exceeds  the  frequency  selective  or  coherence  bandwidth  of  the  chan¬ 
nel.  When  the  scintillations  exhibit  statistical  decorrelation  across  the 
face  of  an  aperture  antenna,  the  channel  may  also  be  referred  to  as  spa¬ 
tially  selective.  Spatially  selective  scintillations  are  therefore  en¬ 
countered  when  the  antenna  aperture  size  exceeds  the  decorrel ation  dis¬ 
tance  of  the  incident  signal. 


Under  conditions  where  the  signal  is  spatially  selective,  the 
antenna  beamwidth  is  smaller  than  the  angle-of-arri val  fluctuations  and 
the  effect  of  the  antenna  is  to  attenuate  the  incident  signal  that  is 
arriving  at  off-boresight  angles.  In  the  spatial  domain,  the  incident 
electric  field  is  somewhat  decorrelated  across  the  face  of  the  antenna. 

The  induced  voltages  in  the  antenna  then  do  not  add  coherently  as  they 
would  for  an  incident  plane  wave  with  a  loss  in  the  gain  of  the  antenna  as 
a  result.  Because  of  this  angular  filtering  or  spatial  selectivity,  the 
second  order  statistics  of  the  signal  at  the  output  of  the  antenna  will  be 
different  than  those  of  the  incident  signal. 

The  effects  of  antennas  on  signals  that  have  propagated  through 
randomly  ionized  media  have  been  reported  by  Wittwer  (1982,  1986)  and 
Knepp  (1983a,  1985).  The  purpose  of  this  report  is  to  review  the  basic 
theory,  starting  with  Maxwell's  equations,  of  radio  frequency  (RF)  signal 
propagation  through  random  media.  Then  several  new  applications  of  an¬ 
tenna  filtering  effects  are  presented.  These  include  the  filtering  of 
anisotropic  signals  with  anisotropic  antennas;  generation  of  realizations 
of  the  impulse  response  function  at  the  output  of  multiple  antennas;  and 
transponder  communication  links  with  two  independent  propagation  paths  and 
four  antennas. 

1.1  THEORY. 

The  starting  point  for  antenna  aperture  effects  calculations  is 
the  generalized  power  spectal  density  (GPSD).  The  first  part  of  this 
report  is  a  review  of  the  derivation  of  the  GPSD.  The  intent  of  this 
review  is  to  give  the  reader  an  understanding  of  the  underlying  physics 
that  are  contained  in  the  GPSD  and  an  understanding  of  the  assumptions 
used  to  calculate  the  GPSD.  The  first  part  of  this  review  follows 
Tatarskii  (1971,  §64-65)  and  the  second  part  follows  Knepp  (1983a). 
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The  derivation  of  the  GPSO  starts  with  Maxwell's  equations  from 
which  the  parabolic  wave  equation  is  derived.  The  parabolic  wave  equation 
can  be  solved  to  give  the  received  electric  field  for  a  specific  electron 
density  distribution  in  the  ionosphere.  However,  the  electron  density 
distribution  is  a  random  process  so  the  received  electric  field  is  also  a 
random  process.  The  parabolic  wave  equation  Is  therefore  used  to  derive 
an  equation  for  the  two-position,  two-frequency,  two-time  mutual  coherence 
function  of  the  electric  field,  r( Ar,Aw,at) .  The  solution  of  the  differen¬ 
tial  equation  for  r,  which  is  a  solution  of  Maxwell's  equations,  then  pro¬ 
vides  a  description  of  the  second  order  statistics  of  the  received  elec¬ 
tric  field.  The  Fourier  transform  of  r  is  the  GPSO  of  the  received 
signal . 

1.2  APPLICATIONS. 

The  second  part  of  this  report  presents  several  new  results. 

The  GPSO  for  anisotropic  scattering  is  used  to  compute  the  mean  power, 
decorrelation  distance,  and  frequency  selective  bandwith  of  the  signal  out 
of  an  anisotropic  antenna.  These  results  allow  an  arbitrary  rotation 
about  the  1 ine-of-sight  between  the  antenna  axis  and  the  natural  axis  of 
the  scattering.  Analytic  results  are  presented  for  antennas  with  Gaussian 
beam  profiles  and  numerical  results  are  given  for  uniformly  weighted 
circular  apertures. 

These  results  are  then  generalized  to  transponder  communication 
link  geometries  where  there  are  two  independent  propagation  paths  through 
disturbed  regions  of  the  ionosphere  and  four  antennas.  It  is  assumed  for 
simplicity  in  these  calculations  that  all  of  the  antennas  have  circularly 
symmetric  antenna  beam  profiles. 


In  the  last  applicat'on  action  analytical /numerical  techniques 
are  described  to  generate  realizations  of  the  impulse  response  function  of 
the  signal  after  propagation  through  randomly  ionized  media  and  reception 
by  multiple  antennas.  These  statistical  realizations  of  the  signal  out  of 
multiple  antennas  have  Rayleigh  amplitude  statistics  and  spatial  and  fre¬ 
quency  correlation  properties  given  by  the  mutual  coherence  function. 

These  realizations  of  the  Impulse  response  function  are  then  used  to  con¬ 
struct  the  received  signal  which  may  be  used  to  exercise  simulations  of 
transionospheric  communications  links. 
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SECTION  2 
THEORY 


The  starting  point  for  antenna  aperture  effects  calculations  Is 
the  generalized  power  spectral  density  (GPSD).  This  section  presents  a 
review  of  the  derivation  of  the  GPSD  and  discusses  the  physics  that  are 
contained  in  this  important  function. 

In  deriving  the  GPSD,  two  key  approximations  are  made  about  the 
spatial  and  temporal  electron  density  fluctuations  that  cause  the  scatter¬ 
ing  in  the  ionosphere.  The  first  of  these  is  the  delta  layer  approxima¬ 
tion  which  says  that  the  scattering  occurs  in  an  infinitesimally  thin 
layer  normal  to  the  1 ine-of-sight.  This  approximation  has  been  relaxed  in 
the  calculations  of  Wittwer  (1982)  and  Knepp  (1983b)  and  has  been  found  to 
result  in  small  errors  in  the  GPSD  provided  that  the  propagation  param¬ 
eters  (frequency  selective  bandwidth,  decorrelation  time  and  decorrelation 
distance)  that  characterize  the  channel  are  properly  specified.  The  delta 
layer  approximation  is  not,  in  general,  adequate  to  calculate  the  propaga¬ 
tion  parameters  (Wittwer  1982). 

The  second  approximation  is  for  the  temporal  variation  model  of 
the  electron  density  fluctuations.  The  usual  approximation  is  Taylor's 
frozen-in  hypothesis  which  treats  the  ionization  fluctuations  or  stria- 
tions  as  rigid  "frozen-in"  structures  that  drift  across  the  line-of-sight. 
Under  this  approximation  there  is  strong  coupling  between  the  spatial  and 
temporal  variations  of  the  random  electric  field  that  is  incident  at  the 
plane  of  the  receiver.  A  second  approximation  proposed  by  Wittwer  (1985), 
called  the  turbulent  approximation,  is  developed  in  this  section  which 
decouples  the  spatial  and  temporal  variations  of  the  random  electric 


fields.  This  approximation  may  be  more  accurate  at  times  before  the  stri- 
ations  have  formed  in  the  ionosphere  or  when  multiple  layers  of  striations 
with  different  drift  velocities  are  in  the  1  ine-of-sight.  The  6PSD  under 
both  of  these  approximations  is  derived  in  this  section. 


An  analytic  solution  is  obtained  in  this  section  for  the  two- 
position,  two-frequency,  two-time  mutual  coherence  function  r(ar,Aai,At)  of 
the  complex  electric  field  incident  on  the  plane  of  the  receiving  antenna. 
This  solution  is  valid  for  arbitrary  1 ine-of-sight  geometries  relative  to 
the  ionization  structures  in  the  ionosphere  that  cause  the  scattering  of 
the  RF  wave.  The  mutual  coherence  function  then  provides  the  basis  for 
the  antenna  aperture  effect  calculations  and  for  the  statistical  signal 
generation  techniques  discussed  in  Sections  3  and  5  of  this  report.  The 
Fourier  transform  of  r  is  the  GPSD  of  the  received  signal. 


2.1  PARABOLIC  WAVE  EQUATION. 

* 

Consider  a  monochromatic  spherical  wave  with  an  electric  field 
<£(r,u),t)  which  is  a  function  of  position  r,  carrier  radi an .  frequency  w, 
arid  time  t.  The  wave  originates  from  a  transmitter  located  at 
r  =  (0,0, -zt)  and  propagates  in  free  space  in  the  positive  z  direction 
until  it  is  incident  on  an  irregularly  ionized  layer  which  extends  from 
0  <  z  <  L  and  is  infinite  in  the  x-y  plane.  After  emerging  from  the  layer 
at  z  =  L,  the  wave  propagates  in  free  space  to  a  receiver  located  at 
r  =  (0,0, z  ).  This  geometry  is  shown  in  Figure  2-1. 


The  propagation  of  the  wave  is  governed  by  Maxwell's  equations: 


$x/+  (1/C)  3J?/3t  »  0  $•€<?  = 

(e/C)  3df/3t  --  0  * 

where  c  s  the  speed  of  light  in  vacuum,  is  the 
the  dielectric  cuiibtant. 


0 

0 

magnetic  field, 


(2-1) 

and  e  is 
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Figure  2-1.  Propagation  geometry. 

The  dielectric  constant  of  the  propagation  medium  undergoes 
random  fluctuations  with  a  characteristic  frequency  which  is  assumed  to  be 
small  when  compared  to  the  carrier  frequency  of  the  wave.  With  this 
assumption,  the  electric  and  magnetic  fields  may  be  written  as  the  product 
of  slowly  varying  complex  envelopes,  denoted  f:  and  H,  times  exp(iu>t): 

/(r.ui.t)  =  !(r  ,(D,t)elu>t 

<V(r,o>,  t)  =  H(r,u),t)eiwt 
Inserting  these  into  Maxwell's  equations  gives 

(2-3a) 
f  2 -3b ) 
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$xf  +  ikH  *  0 
$xh  -  iekf:  =  0 


where  k  =  w/c  is  the  wave  number  of  the  carrier.  After  applying  the  curl 
operator  to  Equation  ?-3a  and  substituting  Equation  2-3b  for  the  $*H 
term,  the  equation  for  ?  becomes 

-  ek2f:  =  0  (2-4) 

The  vector  identity  =  $($•£)  -  v2f  reduces  the  curl  curl  term  in 

Equation  2-4  with  the  result 

$($•?)  -  v2f  -  ek2?  =  0  (2-5) 

The  term  is  reduced  by  expanding  the  divergence  equation  for  £: 

=  0  (2-6) 
or 

=  .  f.^(lne)  .  (2-7) 

The  wave  eauation  for  the  complex  envelope  of  the  electric  field  then 
becomes 

V2?  +  ck2£  +  $[£•$( Ine)]  =  0  (2-8) 

The  dielectric  constant  e  in  a  plasma  at  radio  frequencies  is 
given  approximately  by 

e  =  1  -  u2/u>2  (2-9) 

where  the  plasma  frequency  is 

u>2  =  4n  rac2na('"»t)  •  (2-10) 

p  e  e 
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The  quantity  rg  is  the  classical  electron  radius  (rg=2 .8179*10“15  m)  and 
n  (r ,t)  is  the  free  electron  density  as  a  function  of  position  and  time. 
Equation  2-9  is  valid  when  the  carrier  frequency  is  large  compared  to  the 
plasma  frequency.  The  free  electron  density  is  a  random  variable  that 
will  be  represented  as  a  mean  value  plus  a  random  variation: 

ng(r,t)  =  <r<e>  +  4ng(r,t)  .  (2-11) 

The  electron  density  fluctuation  4ne(r,t)  is  assumed  to  be  a  zero  mean 
random  process  with  a  standard  deviation  that  is  small  compared  to  the 
mean  electron  density.  The  term  ek2  in  the  wave  equation  may  now  be 
rewritten  as 

ek2  =  (2-12) 

where 

]<2  *  (a)/c)2[  l-4nrgc2<ne>/a)2]  (2-13) 

and 

*1  -  [Ane(r,t)/<ne>  ][“>2/(w2-u>p)  ]  .  (2-14) 

The  quantity  Up  is  the  plasma  frequency  evaluated  at  the  mean  electron 
density.  The  magnitude  of  the  gradient  term  $(lne)  in  the  wave  equation 
may  be  estimated  as  follows: 

|$(lnc)|  «  |$[ln(l-«2Ao2)]| 

•  -  (4ure<ng>c2/w2) |$(Ane/<ne>) | 

-  -  ("p/“2)/L0 
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where  L0  is  the  scale  size  of  the  electron  density  fluctuations.  As 
long  as  lQ  »  x,  where  x  is  the  RF  wavelength,  the  term  $[£*^(1ne)]  is 
small  compared  to  ek2!-  and  may  be  ignored.  The  steps  that  follow  are 
therefore  only  valid  when  the  scale  size  of  the  electron  density  fluctua¬ 
tions  is  large  compared  to  the  carrier  wavelength.  With  this  restriction, 
the  wave  equation  becomes 

V2!  +  k*  (l-e^!  =  0  .  (2-15) 

Now  consider  the  complex  components  of  the  electric  field  and 
let 


E(r,w,t)  =  U(r,<n,t)  exp  (-i  /  k  (z')dz')  .  (2-16) 

This  scalar  equation  for  f  may  be  used  because  it  is  usual  for  trans- 
ionospheric  RF  links  to  be  circularly  polarized.  It  is  therefore  not 
necessary  to  carry  out  separate  calculations  for  each  polarization  state. 
The  exponential  term  in  Equation  2-16  represents  the  dispersive  effects  of 
the  smooth  plasma  and  will  be  considered  in  Section  2.9.  The  voltaoe  U 
contains  the  diffractive  effects  that  are  of  interest  under  strong  scat¬ 
tering  conditions.  Substituting  Equation  2-16  for  E(r,u>,t)  in  the  wave 
equation  gives  the  following  differential  equation  for  U: 

vjTJ  +  32U/3z2  -  2i"k 3U/3z  -  k*eiU  *  0  (2-17) 

where 

v|  =  a2/3x2  +  32/3y2  .  (2-18) 

The  complex  amplitude  U  varies  as  the  electron  density  fluctuations.  The 

second  derivative  32 U/ 3 z 2  is  then  the  order  of  U/L2 .  On  the  other  hand, 

the  term  k3U/3z  varies  as  U/XL  .  As  long  as  x  «  L  ,  the  second  derivative 

o  o 

is  small  compared  to  the  first  derivative  and  may  be  iqnored.  This  is 
equivalent  to  neglecting  reflected  rays  and  is  called  the  "parabolic" 
approximation.  The  parabolic  wave  equation  is  then 
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(2-19) 


vju  -  2i k au/Sz  -  k  ElU  =  0 

It  will  be  seen  that  this  parabolic  or  small  angle  scattering  approxima¬ 
tion  is  very  robust  in  that  it  degrades  gracefully  as  the  scattering 
angles  get  large.  The  source  term  exU  in  the  parabolic  wave  equation  is  a 
function  of  the  electron  density  fluctuations  and  of  frequency.  Different 
frequencies  within  a  signal  bandwidth  may  therefore  propagate  along  dif¬ 
ferent  paths  through  the  same  ionization  structure,  '.'.en  this  happens, 
the  propagation  channel  is  said  to  be  frequency  selective. 

2.2  TRANSPORT  EQUATION. 

A  partial  differential  equation  for  the  two-position,  two-fre¬ 
quency,  two-time  mutual  coherence  function  r  is  derived  from  the  parabolic 
wave  eauation  in  this  section  and  is  solved  in  Section  2.5.  This  trans¬ 
port  equation  is  derived  usinq  the  Novikov  theorem  which  requires  that  the 
electron  density  fluctuations  b*3  normally  distributed.  However,  Lee  and 
Jokipii  (1975a)  give  an  alternative  derivation  that  relaxes  this  assump¬ 
tion. 

2.2.1  First  Form  of  the  Transport  Equation. 

The  two-position,  two  frequency,  two-time  mutual  coherence 
function  is  defined  in  a  plane  normal  to  the  1 ine-of-sight  as 

r  =  <U(P1,z,»lft1)U*(p2,2,»2ft2)>  (2-20) 

where  p  is  a  two-dimensional  position  vector  in  the  normal  plane. 

In  order  to  obtain  an  equation  for  r,  the  parabolic  wave  equa¬ 
tion  for  U(pi,z,w1,ti)  is  multiplied  by  U*(P2,z,u'2,t2) .  This  results  in 
the  following  equation: 
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(l/kl)^12U(Pl,ZIa)l,t1)U*(P2,Z,u>2,t2)  - 

2i [3U(Pi,Z,Ui,ti)/3z]U*(P2.2.“2tt2)  - 

kici(Pi.z,<ui,t1)U(Pifz,u1,t1)U*(p2,z,“2»t2)  s  0  (2-21) 

where  kj  is  given  by  Equation  2-13  evaluated  at  Wj,  e(pj »2»uj »tj )  is  given 
by  Equation  2-14,  and  V.2  is  given  by  Equation  2-18  evaluated  at  p^_.  A 
similar  equation  can  be  written  down  by  interchanging  the  subscripts  1  and 
2  and  by  taking  the  complex  conjugate  with  the  result: 

(l/k2)^2U*(P2,Z,u)2,t2)U(Pi,Z,w1,t1)  + 


2i  [3U*(p2,Z,a)2,t2)/3z]U(p1,Z,u1,t1)  - 

k2el(02,Z,w2,t2)U*(P2,Z,a>2,t2)U(pi,Z,u)1,t1)  ~  0  •  (2-22) 

Upon  subtracting  Equation  2-22  from  Equation  2-21  and  taking  the 
expectation  value,  the  equation  for  r  is 

(l/ki)^2r  -  (l/k2)^_L22r  -  2i ar/az  - 

ki<ei(pi,z,w1,ti)U(Pi,z,^i,ti)U*(p2,z,o>2,t2)>  + 
k2<ei(P2»z,u,2.t2)U(Pi,z,u>i,ti)U*(P2,z,w2,t2)>  =  0  .  (2-23) 


The  expectation  of  the  two  source  terms  in  this  equation  must  be 
carefully  evaluated.  They  involve  the  product  of  UU*  and  where  t\  is 
proportional  to  the  fluctuations  in  the  electron  density.  However,  the 
electric  field  complex  envelope  U  is  a  function  of  the  electron  density 
fluctuations  that  are  encountered  along  the  propaaation  path.  Therefore  U 
and  gi  are  correlated. 


*  •  -  ’s  • v  M  v  j. *  •  - 
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2.2.2 


Novikov  Theorem. 


The  Novikov  theorem  is  used  to  evaluate  the  source  terms  in 
Equation  2-23.  This  theorem  is  proved  in  Tatarskii  (1971,  §65)  and 
Ishimaru  (1978,  pp.  457-458).  The  theorem  states  that 

<f(h*[f)>  a  /  <f ($)f (ft*  )><4$[f  ]/6f ($'  )>  dn  &'  (2-24) 

where  f($)  is  a  zero-mean,  normally  distributed  random  function  of  the  n 
dimensional  vector  $  and  S^/Sf  is  a  functional  derivative.  In  applying 
this  theorem,  f  =  and  ♦  =  UU*.  The  theorem  is  proved  by  expanding  <t>(f) 

in  a  Taylor  series. 

2.2.3  Source  Terms. 

Before  proceeding  with  evaluatina  the  source  terms,  it  will  be 
convenient  to  write  as  the  product  of  a  frequency  term  and  a  term  that 
varies  only  with  space  and  time: 

ei(p,z,w,t)  «  5(p,Z,t)B(u>)  (2-25) 

where 

£(p,z,t)  =  Ane(p,z,t)/<ne>  (2-26) 

is  a  random  function  of  the  electron  density  fluctuations  and 

8(a))  =  ^/(u2-®*)  (2-27) 

is  a  deterministic  function  of  frequency  and  the  mean  free  electron 
dens i ty. 


A  straightforward  application  of  the  Novikov  theorem  qives  for 
the  first  source  term  in  Equation  2-23 
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Si  =  k i S  (o>i  )<C(pi  ,z ,  t 1  )U(  Pi  ,z,wi  ,t  i  )U*(p2  ,z,<*>2  >^2  )> 


=  MNi)  /  dz'  //  d2o'  /  dt '  <C(Pi,2,t1)5(p',z,,t,)>  x 
.00  *  00 

<[au(pi,z,0)1,ti)/65(pl  ,z*  ,t,)]U*(p2,Z,W2»t2)  + 

U(Pi,Z,0)i,t1)[6U*(p2,z,a)2,t2)6C(?,,z,,f)]>  .  (2-28) 

At  this  point,  the  electron  density  fluctuations  are  assumed  to 
be  stationary  and  delta-correlated  along  the  z  axis.  This  Markov  assump¬ 
tion  has  the  mathematical  form 

<C(P,z,t)5(p',z' ,t')>  =  «F(z-z' )A(p-p* , t-t ’ )  .  (2-29) 

The  notation  5p(x)  is  used  throughout  this  report  for  the  Dirac  delta 
function  to  distinguish  6p(x)  from  the  functional  differential  <5<i>/<5f  and 
from  the  parameter  5  defined  in  Section  2.6.  The  Markov  assumption  is 
discussed  in  some  detail  by  Tatarskii  (1971,  §64)  and  is  based  on  the  fact 
that  fluctuations  in  the  dielectric  constant  in  the  direction  of  propaga¬ 
tion  have  little  effect  on  the  transverse  fluctuation  characteristics  of 
the  electric  field.  It  is  the  fluctuations  of  the  dielectric  constant 
transverse  to  the  direction  of  propagation  that  dominate  the  scattering 
and  the  transverse  fluctuations  of  the  electric  field. 

Under  the  assumption  of  small  angle  scattering  for  which  the 
parabolic  wave  equation  is  valid,  the  component  of  the  electric  field 
traveling  in  the  backward  direction  will  be  negligible  compared  to  the 
component  of  the  electric  field  traveling  in  the  forward  direction.  The 
electric  field  U(p,z,<u,t)  may  then  be  assumed  to  depend  on  5(p',z',t') 
only  for  z'  <  z  (i.e.  the  electric  field  does  not  depend  on  electron  den¬ 
sity  irregularities  that  have  not  yet  been  encountered  along  the  forward 
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propagation  path).  Also,  U(p,z,«i,t)  depends  on  C(p',z',t')  only  for 
t'  <  t  (i.e.  the  electric  field  does  not  depend  on  irregularities  that 
have  not  yet  occurred).  Thus  6U(p,z,u),t)/65(p' ,z' ,t' )  =  0  for  z'  >  z  and 
for  t'  >  t. 

The  source  function  Si  may  now  be  rewritten  as 
z  *  t 

Si  =  k i6 1  /  dz'  6 p( z-z 1 )  //  d2p'  /  dt '  A( pi-p ' , 1 1 -t ' )  * 

-  ao  wa>  «*oo 

<[«U(Pi,Z,a»i,ti)/«C(p'  , Z 1  ,t')  ]U*(p2,Z,u>2,t2)  + 

U(Pl,Z,<*)Uti)[6U*(p2,Z,ai2,t2)/65(p,,Z,,t,)]>  (2-30) 

where  S  .  is  given  by  Equation  2-27  evaluated  at  w . .  Recalling  that 

J  J 

a 

/  5  (x-a)dx  =  1/2  (2-31) 

^00 

further  reduces  the  source  term  Si  to 

«  t 

Si  =  ( k 1 0 1/2 )  //  d2P'  /  dt'  A(Pi-P,,ti-t')  x 

.00  .  oo 

<[6U(Pi,Z,wi,ti)/55(p,,Z,t’) ]U*(p2,Z,u>2,t2)  + 

U(Pi,Z,w1,ti)[6U*(p2,z,w2,t2)/'55(p,,z,t')]>  .  (2-32) 
Tne  SU/SC  terms  in  this  expression  will  be  evaluated  next. 

The  parabolic  wave  equation  is  used  to  evaluate  the  functional 
derivatives  6U/6£.  Integrating  this  equation  from  -•  to  z  results  in 


m 


m 


M 


m 


$ 


/  ^  (j(p,c,u,t)dc  -  2ik[U(p,z,u,t)  -  UQ{ P ,io) ]  - 


k2e(u>)  j  €(pt2,t)U(p,2,«tt)dc  =  0 


(2-33) 


where  U0(p,w)  is  the  transmitted  signal.  After  applying  the  operator 
S/$£(p' ,2' ,t' )»  where  -®  <  z1  <  z  and  -®  <  t'  <  t,  and  noting  that 

«€(p,z.t)/«5(p,.2*  .f)  *  6F(z-z')«F(5-P,)«F(t-f)  ,  (2-34) 

Equation  2-33  becomes 

2  ik  5U(p,Z,i*»ft)/65(p'  ,z'  ,t' )  +  k  6(tii)5fr(  p-p ' )  6p(  t-t ' )  U(  P,Z '  ,iu,t )  + 

2  _2 

/  [k  6(w)c(p,c,t)  -  ^|][5U(p,;,a.,t)/65(;',z,,f)]dc  =  0  .  (2-35) 

z'  ± 

The  lower  limit  of  the  integral  in  this  equation  is  z'  because 
<5U(p,z,o),t)/6C(p,,zl  ,t 1 )  is  zero  for  z  <  z'.  Because  the  source  term  Si 
contains  factors  of  the  form  SU(p,z,w,t)/45(p',z,t') ,  z'  is  set  equal  to 
z  in  Equation  2-35  with  the  result 


<5U(p,Z,n>,t)/S5(p' ,z,t')  =  (  ik8(a>)/2)'Sp(p-p,)6F(t-t,)U(p,z,u,t)  (2-36) 


A  similar  expression  may  be  written  down  for  the  second  source  term  in 
Equation  2-23: 

=  (iki8182/4)A(Pi-p2,t1-t2)  r  -  ( ik2822/4)A(($,0)  r  .  (2-38) 

2.2.4  Second  Form  of  Transport  Equation. 

The  transport  equation  for  the  mutual  coherence  function  is  now 
given  by  combining  Equations  2-37  and  2-28  with  2-23  with  the  result 

ar/r2  +  (i/2)[(l/kl)v^  -  (i/k2)vj]  r  - 

(l/8)[2k1k26l82A(p1-P2.t1-t2)  -  (k12612+k2262Z)Ao]  T  (2-39) 
where  Aq  =  A(5,0). 

The  differential  equation  for  T  will  be  solved  by  first  lettinq 
r  =  r  r:  Where  r  is  the  free  space  solution.  The  well  known  solution  of 
the  wave  equation  for  the  electric  field  in  free  space  (Equation  2-15  with 
ei  -  0)  may  be  written  down  directly.  The  Fresnel  approximation  that 
z  »  |?|  is  then  used  to  expand  the  electric  field  and  the  free  space 
solution  rQ  is  computed.  The  quantity  Tq  contains  the  1/z2  term  that 
partly  determines  the  mean  power  at  the  receiver.  The  next  step  is  to 
derive  a  differential  equation  for  from  Equation  2-39  and  the  free 
space  solution.  It  is  the  mutual  coherence  function  Tj  that  determines 
the  second  order  statistics  of  the  received  signal. 

2. 2. 4.1  Free  Space  Solution  r0.  in  free  space  and  for  spherical 
geometry,  the  complex  envelope  of  the  electric  field  is 


It  is  easy  to  verify  that  this  is  a  solution  of  the  wave  equation  with  e-i 
set  to  zero.  Under  the  assumption  of  small  angle  scattering,  z2  »  x2  +  y2 
in  the  region  of  Interest  and  |r|  may  be  expanded  as 

j r |  «  (x2+y2+z2) 1/2  -  z  +  (x2+y2)/2z  .  (2-41) 

With  this  Fresnel  approximation, 

?  *  (?q/z)  exp  [-1kz-1k(x2+y2)/2z]  .  (2-42) 

After  recalling  that  U  =  exp(ikz)!,  the  free  space  mutual  coherence 
function  is 

rQ  =  <U(Pi,Z,«i)U*(P2,Z,“2)> 

=  z*2  exp  [ ik i (x i2+y i2 )/2z  +  ik2(x22+y22 )/?z]  •  (2-43) 

2. 2. 4. 2  Differential  Equation  for  r1#  After  substituting  r  *  rQr1  into 
Equation  2-39  and  using  Equation  2-43,  the  equation  for  becomes 

3lV3z  »  (i/2k1)7i2r1  -  ( i /2k2) VjJrj  + 

(xi/z)ar1/3xl  +  (yi/z) ar1/3y1  +  (x2/z)3r!/3x2  +  (y2/z)8ri/3y2 
-  s  rx  =  0  (2-44) 

where  the  source  term  is 

S  =  [2k1k28l62A(p1-p2,t i-t2)  -  (k 12S12+k22622 )A0]/8  •  <2'45) 

In  order  for  Tj  to  represent  a  statistically  stationary  random 
process  in  space,  frequency,  and  time,  Fi  must  be  a  function  only  of  the 
differences  Pi  -  p2»  *  ^z*  snd  tj  -  t.2.  It  is  therefore  useful  to 

transform  Equation  2-44  to  sum  and  difference  spatial  and  frequency 
coordinates 
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X  s  ( x  i  +x  2 )  /  2 


y  «  (yi+y2)/2 
5  ‘  X!  ■  x2 

n  =  yi  -  y2 

72  ,  a2/3x2  +  32/3y2  (2-46) 

v2  *  32/3;2  +  32/3n2 
d 

$  =  3 2 / 3 X 3 C  +  32/3Y3ti 

s  d 

k$  s  ( k i  +k 2 )/2 

^d=  k 1  ”  ^2 

After  some  manipulations,  the  equation  for  r2  reduces  to 

3IV3Z  -  (i/2)(k2-k2/4)"1  (kdV2  +  (kd/4)72  -  k^s^d]  Tj  + 

Z_1[X  3/3X  +  Y  3 / 3 Y  +  C  3/3C  +  n  3/»n]  -  SI^  =  0  .  (2-47) 

The  boundary  condition  on  is  that  r^z^-z^)  =  1  independent 
of  the  other  spatial  coordinates.  Also,  the  source  term  S  under  most  con¬ 
ditions  is  a  function  only  of  the  difference  coordinates.  It  will  there¬ 
fore  be  assumed  that  r2  Is  Independent  of  X  and  Y.  However,  the  source 
term  will  be  a  function  of  X  and  Y  if  the  spatial  extent  of  the  scattering 
region  is  small  as,  for  example,  in  a  barium  cloud.  The  assumption  that 
r j  is  independent  of  X  and  Y  then  requires  that  the  disturbed  region  in 
the  Ionosphere  be  large  compared  to  the  region  from  which  scattered  signal 
energy  is  received. 
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The  transport  equation  may  be  further  reduced  by  noting  that  the 
1  h  term,  when  z  is  large,  will  be  small  compared  to  the  other  terms  and 
this  term  may  be  neglected.  The  differential  equation  for  rx  then  becomes 

3IV3Z  -  (ik(j/2)(k*-k^/4)~1v*r1  -  srx  -  0  .  (2-48) 


2.3  OELTA  LAYER  APPROXIMATION. 

As  an  RF  wave  propagates  through  a  thick,  irregularly  structured 
ionization  layer,  the  wave  first  suffers  random  phase  perturbations  due  to 
random  variations  in  the  index  of  refraction.  As  the  wave  propagates 
farther,  diffractive  effects  introduce  fluctuations  In  amplitude  as  well 
as  phase.  If  the  standard  deviation  of  the  phase  fluctuations  that 
are  suffered  by  the  wave  is  large,  then  the  amplitude  fluctuations  are 
characterized  by  a  Rayleigh  probability  distribution  when  the  wave  emerges 
from  the  layer.  The  delta  layer  approximation  assumes  that  the  phase  and 
amplitude  fluctuations  are  imparted  on  the  wave  in  an  infinitesimally  thin 
layer.  This  assumption  is  consistent  with  the  Markov  assumption  made  in 
deriving  the  differential  equation  for  the  mutual  coherence  function. 

An  analytic  solution  for  the  two-position,  two-frequency  mutual 
coherence  function  has  been  obtained  by  Wlttwer  (1979)  and  Knepp  (1983b) 
for  a  thick  ionization  layer.  The  analytic  form  of  this  solution  is  suf¬ 
ficiently  complex,  however,  that  the  necessary  Fourier  transforms  required 
to  compute  the  GPSD  cannot  be  performed  in  closed  form.  The  complex  ana¬ 
lytic  form  is  simplified  by  the  use  of  the  delta  layer  approximation  to 
obtain  tractable  expressions  for  the  mutual  coherence  function  and  the 
GPSD.  Wittwer  and  Knepp  have  evaluated  the  accuracy  of  the  delta  layer 
approximation  in  detail  as  it  affects  the  delay  distribution  of  the  re¬ 
ceived  signal  and  have  found  that  the  maximum  error  Is  small  for  trans- 
ionospheric  satellite  communication  link  geometries  as  long  as  the  param¬ 
eters  of  the  GPSO  are  properly  selected.  Wittwer  (1979)  has  derived 
expressions  for  the  signal  parameters  of  the  GPSD  that  include  the  effects 
of  a  thick  scattering  layer. 
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Now  a  relationship  between  the  electron  density  fluctuations  and 
the  phase  variations  imposed  on  the  wave  may  be  calculated.  The  differen¬ 
tial  phase  change  of  the  wave  along  the  propagation  path  i  is 

dVdi  =  reXAne(p,z,t)  •  (2-49) 

Under  the  assumption  of  small  angle  scattering  used  to  derive  the  differ¬ 
ential  equation  for  r,  d/d!  -  d/dz  and  the  total  phase  change  of  the  wave 
is 

=  rgX  /  Ane(p,z,t)dz  =  rex<ne>  /  C(p,z,t)dz  ,  (2-50) 

Integrated  through  the  Ionization  layer.  The  autocorrelation  function  of 
the  phase  *  is 

<<t>(p,t)4>(p',t'  )>  =  (rex<ne>) 2  /  dz  /  dz’  <5(p,z,t)5(p'  ,z' , t '  )> 
s  (rex<ne>) 2  /  dz  A(p-p'.t-t') 

=  (reX<ne>)  2  L5A(p-p',t-f)  (2-51) 

where  L<s  is  the  delta  layer  thickness.  The  Markov  approximation 
(Equation  2-29)  has  been  used  in  evaluating  the  autocorrelation  of  5. 
However,  it  is  shown  in  Appendix  A  that  Equation  2-51  is  valid  as  long  as 
the  scattering  layer  thickness  is  large  compared  to  the  parallel  decor¬ 
relation  distance  of  the  electron  density  fluctuations.  Finally,  the 
phase  variance  imparted  on  the  wave  is 

•  <'e1<V)  Vo  •  (2-52) 

The  quantity  A0  depends  on  the  power  spectrum  of  the  electron  density 
fluctuations  in  the  ionosphere.  The  value  of  A0  for  a  three-dimensional 
k-4  in  situ  power  spectrum  and  for  the  delta  layer  approximation  is  given 
in  Appendix  B. 


21 


In  general,  only  part  of  the  total  phase  variance  results  in  the 
Rayleigh  amplitude  fluctuations  that  are  of  described  by  the  GPSD. 

Wittwer  (1979,  1980)  calls  this  part  the  Rayleigh  phase  variance.  The 
rest  of  the  total  phase  variance  is  associated  with  the  mean  dispersive 
effects  described  by  the  exponential  term  of  Equation  2-16.  Physically, 
it  is  the  smaller  sized  electron  density  fluctuations  that  result  in  dif¬ 
fractive  effects  and  the  larger  sized  fluctuations  that  result  in  disper¬ 
sive  effects.  Both  Wittwer  and  Knepp  (1982)  describe  how  the  Rayleigh 
phase  variance  may  be  separated  from  the  total  phase  variance.  In  the 
developments  that  follow,  the  phase  variance  in  Equation  2-52  will  be 
assumed  to  be  the  Rayleigh  phase  variance  associated  with  amplitude 
fluctuations. 

2.4  FROZEN-IN  APPROXIMATION. 

Under  the  frozen-in  approximation,  the  temporal  variation  of  the 
electron  density  fluctuations  is  qiven  by 

ne(P,2,t)  *  ng(p-vt,z,0)  .  (2-53) 

This  equation  is  valid  if  the  electron  density  fluctuations  with  a  scale 
size  L0  do  not  appreciably  change  their  shape  within  the  time  required 
for  the  structures  to  drift  a  distance  L0.  This  is  called  Taylor's 
frozen-in  approximation.  With  this  approximation,  the  function  A  becomes 

A(p,t)  =  A( p- vt ,0)  .  (2-54) 

This  approximation  is  accurate  for  ionospheric  conditions  where 
the  ionization  has  broken  up  into  a  single  layer  of  striations  aligned 
with  the  earth's  magnetic  field  lines.  The  frozen-in  approximation  may 
not  be  valid  before  striations  have  formed  or  when  there  are  multiple 
scattering  layers  in  the  1 ine-of-sight.  Alternative  forms  for  the 
coherence  functionand  the  generalized  power  spectral  density  for  this 
latter  situation  are  considered  in  Section  2.9. 
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2.5 


SOLUTION  OF  THE  TRANSPORT  EQUATION. 


Before  proceeding  with  the  solution  of  the  transport  equation, 
it  is  convenient  to  expand  the  source  term  S  by  making  two  non-restricti ve 
assumptions.  First,  it  will  be  assumed  that  the  RF  frequencies  of 
interest  are  large  compared  to  the  plasma  frequency.  Second,  it  will  be 
assumed  that  k<j  is  much  smaller  than  ks.  The  solution  obtained  will 
then  be  valid  for  a  small  range  of  frequencies  around  ks  and  for  fre¬ 
quencies  large  compared  to  typical  peak  ionospheric  plasma  frequencies  of 
a  few  hundred  MHz.  With  these  assumptions,  the  source  term  becomes 

S  •  (k“/4k^)[A(Sd,td)  -  A0]  -  (k^/8k‘)Ao  (2-55) 

where  is  a  two  dimensional  position  vector  in  the  t-n  plane,  t^  is  the 

time  difference,  and  k  =  u  /c. 

P  P 

The  last  term  in  Eauation  2-55  is  a  function  of  freauency  and  z 
but  is  independent  of  t  and  n.  This  suggests  that  another  useful  factor¬ 
ization  is  l*!  =  r2r3  where  r3  is  independent  of  C  and  n.  After  substi¬ 
tuting  r2r3  for  in  Equation  2-48  and  separating  variables,  the  result 
is 

(l/r2)ar2/3z  -  ( ik d/2k |r2 ) v2 r2  -  (k^/4kj) [A(pd,td)  -  aq]  + 

(l/r3)3r3/3Z  +  (kd2kj/8kj)  Aq  =  0  .  (2-56) 

The  last  two  terms  of  this  equation  are  spatial  functions  only  of  z  and 
therefore  must  be  separately  equal  to  zero  for  arbitrary  values  of  c  and 
n. 


The  source  term  in  the  r3  equation  is  only  non-zero  within  the 
delta  layer.  Thus  from  the  transmitter  to  the  delta  layer,  r3  is  unity. 
With  this  boundary  condition,  the  solution  of  the  r3  equation  is 


(2-57) 
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r3  =  exp  [-(kpk^/8kJ)AQ(z-zt)]  ,  z  >  zt  . 

This  term  gives  the  effect  of  the  different  transit  times  of  the  different 
frequencies  due  to  the  frequency  dependence  of  the  index  of  refraction. 

Now  the  equation  for  r2  may  be  solved  using  the  delta  layer 
approximation.  The  equation  for  r2  is 

(l/r2)3r2/3z  -  (ikd/2k2r2)?2r2  -  (kJ/4k2)[A(pd,td)  -  aq]  «  0  .  (2-58) 

Within  the  delta  layer,  the  k<j  term  is  small  compared  to  the  source  term 
and  may  be  ignored.  This  is  equivalent  to  saying  that  diffractive  effects 
are  not  important  within  the  delta  layer  (Lee  and  Jokipii  1975b). 
Integrating  Equation  2-58  through  the  delta  layer  gives  the  value  of  r2  at 
the  point  where  the  wave  emerges  from  the  delta  layer: 


r2  =  exp  i(kpL6/4ks2)[A(pd,td)  -  Ao]}  .  (2-59) 

The  solution  from  this  point  proceeds  as  follows.  Between  the 
delta  layer  and  the  receiver,  the  signal  propagates  in  free  space  so  the 


equation  for  is 

arj/az  -  ( ik(j/2k2)v2r1  =  o  (2-60) 

The  v^!  term  in  this  equation  gives  the  effects  of  diffraction  on  the 
signal  as  it  propagates  from  the  delta  layer  to  the  receiver.  The 
boundary  condition  is  =  r2r3  at  the  point  where  the  wave  emerges  from 
the  delta  layer.  This  equation  is  easily  solved  by  taking  the  Fourier 
transform  from  spatial  coordinates  Z  and  n  to  anqular  coordinates.  First 
it  is  convenient  to  transform  variables  to 

u  =  z/z 

v  =  n/z 
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The  angular  variables  and  are  then  independent  of  2.  After  the 
change  in  variables,  the  equation  for  1^  becomes 


31^/az  -  (ikd/2k2z2)(32/3u2  +  32/ 3v2 )  r 1  *  0  .  .  (2-61) 


The  Fourier  transform  pairs;  from  spatial  coordinates  p  in  the 
plane  normal  to  the  1 ine-of-slght  to  angular  coordinates  from  carrier 
frequency  differences  u>d  to  delay  t,  and  from  time  differences  t^  to 
Doppler  frequencies  u^;  are  defined  in  this  report  to  be 


r(i^)  =  //  exp  (-iit^*p)r(p)d2p 


(2-62a) 


r(t)  =  (2tt)*1  /  exp  (ioidT)r(iuc1)dujd 


(2-62b) 


and 


r(a)Q)  =  /  exp  (-iu>0td)r{td)dtd 


r(p)  =  (2it)-2  //  exp  (i^*p)r(^)d2^ 


r(^d)  =  /  exp  (-iadT)r(u»d)dT 


(2-62c) 


( 2-63a) 


( 2-6  3b ) 


s 

I 

! 

a 


X5 


is 


r(td)  *  (2it )_1  /  exp  (i«0td)r(td)d«D 


( 2-6  3c ) 


Upon  transforming  from  u  and  v  to  Ku  and  K v,  the  equation  for  Tj 


ar , /3Z  +  (ik  /2k2z2)(K2+K2)r,  =  0 
d  s  u  v 


(2-64) 
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This  equation  is  integrated  from  z  «  zt  +  to  z  *  +  zr  which  gives 

Fl  (Ku’Kv’Zt+zr,aid,td)  =  r'1^Ku’Kv’2t+L6,a,d,td^  x 
exp  [-(ikd/2k2)(K2+Kj)Y]  (2-65) 

where 


Y 


zt+zr 
I  z-2dz 

zt+t-6 


(zt+zr)(zt+Lfi) 


(2-66) 


The  expression  for  y  may  be  simplified  by  setting  L5  to  zero  which  is 
consistent  with  the  delta  layer  approximation. 


2.6  QUADRATIC  PHASE  STRUCTURE  APPROXIMATION. 


At  this  point,  a  formal  solution  of  the  Fourier  transform  of  the 
two-Dosition,  two-frequency,  two-time  mutual  coherence  function  at  the 
receiver  has  been  computed  in  terms  of  the  structure  function  MPd,td)  - 
A0.  As  a  practical  matter,  the  quadratic  phase  structure  approximation 
is  reauired  to  make  the  exponent  in  Equation  2-59  quadratic  In  the  spatial 
variables.  This  gives  the  resulting  coherence  function  a  tractable  mathe¬ 
matical  form.  However,  for  the  small  angle,  strong  scattering  conditions 
considered  in  this  report,  the  correlation  distance  of  the  signal  will  be 
much  smaller  than  the  correlation  distance  of  the  electron  density  fluc¬ 
tuations.  The  coherence  function  will  then  be  determined  primarily  by  the 
values  of  A(Pd, td)  -  Aq  at  small  values  of  Pd  and  td.  A  Taylor  series 
expansion  of  A(pd,td)  -  A  ,  keeping  only  the  quadratic  terms,  will  then 
provide  a  reasonably  accurate  mutual  coherence  function. 


Under  the  frozen-in  approximation,  the  temporal  variation  of 
A( pd ,t<j )  is  9^ven  by  A(p^-vtd,0).  It  is  therefore  necessary  to  consider 
only  the  spatial  form  A( Pd ) .  The  Taylor  series  expansion  will  make  the 
following  calculations  independent  of  the  functional  form  of  A(pd)  as  long 
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as  the  second  derivative  of  A  exists  for  pd  equal  to  zero.  A  detailed 
discussion  of  A{pd)  will  therefore  be  deferred  to  Appendix  B.  However, 
some  generic  properties  of  A{pd)  need  to  be  considered  here  in  order  to 
specify  the  functional  form  of  the  quadratic  terms  in  the  Taylor  series 
expansion. 

The  propagation  coordinate  systems  are  shown  in  Figure  2-2.  The 
z  axis  is  along  the  1 ine-of-sight  and  the  t  axis  is  aligned  with  the  geo¬ 
magnetic  field  lines  §  at  the  elevation  of  the  delta  layer.  The  r  and  s 
axes  are  in  a  plane  normal  to  the  magnetic  field.  The  penetration  angle  \|» 
is  the  angle  between  the  z  axis  and  the  §  axis.  In  the  plane  of  the 
receiver,  the  x  axis  direction  is  given  by  the  cross  product  of  the  $ 

A 

vector  and  the  z  unit  vector. 

Consider  now  the  functional  form  the  power  spectrum  of  the 
electron  density  fluctuations  in  the  r-s-t  coordinate  system.  The  power 


x  =  { x  z/sirvf» 


Figure  2-2.  Propagation  coordinate  systems. 


spectrum  is  usually  assumed  to  be  a  function  of  the  quantity  L2K2  +  L2K2 
L2K*  where  lr,  L$ ,  and  are  the  scale  sizes  of  the  fluctuations  in  the 
three  directions.  The  power  spectrum  must  be  rotated  into  the  x-y-z 
coordinate  system  in  order  to  calculate  A(p^).  Before  performing  this 
rotation,  the  usual  assumption  that  the  electron  density  fluctuations  are 
elongated  in  the  t  direction  and  are  symmetric  about  the  t  direction  will 
be  made.  The  scale  sizes  are  then 


Lr  *  Ls  =  Lo 
Lt  *  ^Lo 


(2-67) 


where  q  (q  2  1)  is  the  axial  ratio.  After  the  rotation,  the  power 
spectrum  will  be  a  function  of  the  quantity  (Wittwer  1979)  l*K2  +  + 

LzKz  +  2LyzK/z  Where  7 

Lx  =  Lo 


=  Lo(cos2'l'+q2sin24’) 1/2  =  Lq/6 

Lz  =  L0(sin2^+92cos2'i') 1//2 
Lyz  =  L0(9?-l)sin4-cos'l' 


(2-68) 


The  parameter  6  (1/q  £  6  £  1)  appears  throughout  subsequent  sections  of 
this  report  and  is  defined  to  be  the  ratio  of  the  x  and  y  scale  sizes: 


5  =  l^/L  =  (cos  2<P+q  2sin24>)~  1//2 


(2-69) 


/,  the  function  A(p^)  is  the  two  dimension  Fourier  transform 


of  the  power  spectrum  4^: 


Setting  Kz  5  0  in  the  power  spectrum  in  this  equation  results  in  having 

be  a  function  of  L2K2  +  L2K2  and  in  having  A(p  .)  calculated  in  the 
s  a  x  j  y  u 

z  s  0  plane  of  the  delta  layer.  After  the  Fourier  transform  is  performed, 
A  will  be  a  function  of  the  quantity  x2/L2  +  y2/l2  (see  Appendix  B)  or 
equivalently  of  (x2+<52y2)/L2 . 


The  quadratic  phase  structure  approximation  now  takes  the  form 

A(pd)  *  A0[l-A2(x2+62y2)/L2]  .  (2-71) 

The  coefficient  A2  is  calculated  in  Appendix  B  for  a  K‘4  electron  density 
fluctuation  spectrum.  The  coefficient  A0  is  given  in  terms  of  the  phase 
variance  o2  by  Equation  2-52. 
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2.7  MUTUAL  COHERENCE  FUNCTION. 

The  mutual  coherence  function  r.j  at  the  plane  of  the  receiver  is 
computed  as  outlined  in  Section  2.5.  At  this  point,  the  frozen-in  approx¬ 
imation  will  be  used  to  include  explicit  temporal  variations.  In  order  to 
do  this,  the  drift  velocity  must  be  specified.  To  be  conservative,  the 
drift  velocity  will  be  assumed  to  be  along  the  direction  with  the  smallest 
scale  size  for  the  electron  density  fluctuations  (i.e.  the  x  direction). 

The  boundary  condition  for  rx  at  the  delta  layer  is  aqain  rz  = 
r2r3  evaluated  at  z  =  l6.  Using  Equation  2-52  to  relate  Aq  to  the  phase 
variance  and  then  setting  the  delta  layer  thickness  L<$  to  zero  gives  the 
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following  for  r2  and  r3  at  the  point  where  the  wave  emerges  from  the  delta 
layer.  From  Equation  2-59  using  the  quadratic  phase  structure 
approximation,  the  diffraction  term  boundary  condition  is 

r,  =  exp  {-o2A2[U-vxtd)2  +  52n2]/L2}  (2-72) 
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(2-73) 


and  from  Equation  2-57,  the  retraction  term  boundary  condition  is 
r3  =  exp  [-°J«*>J/(2a>2)] 

where  u)q  =  ck$  is  the  carrier  radian  frequency.  Upon  changing  variables 
in  the  r2  equation  to  the  dimensionless  spatial  coordinates  u  and  v  and 
performing  the  Fourier  transform  to  and  Kv  coordinates,  the  boundary 
condition  for  is 

MKu’WvV  =  [*LS/(a|A2zt6)]  X 

exp  [- /(2«q) ]  exp  [ 1Ku<vxtd/zt ) ]  x 

exp  [-L2(K^K2/62)/(4a|A22|)  ]  .  (2-74) 

The  solution  at  the  plane  of  the  receiver  is  given  by  Equation  2-65  and 
this  boundary  condition. 

Aftjr  performing  the  inverse  Fourier  transform  on  the  solution 
ri(Ku,Kv,t'LJzr,u)d,td)  and  converting  to  unnormalized  distance  units  x  and 
y,  the  mutual  coherence  function  may  be  written  as 

Fl(x,y,o>d,td)  *  [(l+iAa)d/aicoh)(l+iA62t0d/(0coh)]-1/2  * 
exp  [-a2a>2/(2w2)]  x 
exp  [-(x-vetd)2^-2/(l+iAa)d/oJcoh)]  x 

exp  [-62y2Z-2/(l+iA62wd/u>coh)]  (2-75) 

where 

A  =  [2/U+64)]1/2  (2-76) 

and 

ve  =  \<V2r>/2t  <2'77> 

is  the  effective  drift  velocity  of  the  ionosphere  as  seen  at  the  receiver. 
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The  decorrelation  distance  £0  of  the  received  electric  field 
is  defined  as 

£o  =  { Vzr)Lo/(zt/A2  V  (2_78) 

and  the  coherence  bandwidth  <»COh  is  defined  as 

“coh  *  *“JLo{2t+,r)/(2ccJA»ztlr)  •  <2-79> 

It  will  be  shown  below  that  iQ  agrees  with  the  usual  definition  of  de- 

correlation  distance  and  it  will  be  seen  in  the  next  section  that  the 

coherence  bandwidth  is  proportion  to  2*f0  where  f  is  the  frequency 

selective  bandwidth  of  the  signal.  It  is  clear  from  the  form  of  Tj  in 

Equation  2-75  that  the  coherence  bandwidth  could  have  been  defined  as 

“coi/-  The  ^actor  A  has  been  deluded  here  in  order  to  simplify  the 

relationship  between  w  and  f  . 

coh  o 

Equations  2-78  and  2-79  are  only  valid  under  the  delta  layer 
approximation  and  do  not  reflect  now  these  parameters  are  actually  calcu¬ 
lated.  Wittwer  (1979,  1980)  has  derived  expressions  for  the  decorrelation 
distance  and  the  coherence  bandwidth  that  are  valid  for  more  general  scat¬ 
tering  layer  geometries.  It  is  the  expressions  of  Wittwer  that  are  used 
in  signal  specifications  to  calculate  these  signal  parameters.  This  is 
also  true  for  the  decorrelation  time  which  is  formally  defined  under  the 
delta  layer  and  frozen-in  approximations  in  the  next  section. 

For  satellite  communications  links,  the  distance  from  the  iono¬ 
sphere  to  the  satellite  is  typically  much  larqer  than  the  distance  from 
the  ionosphere  to  the  around.  If  the  transmitter  is  on  the  satellite  then 
zt  >  zr,  and  if  the  transmitter  is  on  the  ground  then  zt  <  zp.  8ecause 
the  expression  for  the  decorrelation  distance  is  not  reciprocal  in  zt  and 
zr,  the  value  of  lQ  depends  on  the  direction  of  propagation.  However,  the 

expression  for  the  coherence  bandwidth  is  reciprocal  in  z,  and  z  so  w  , 

r  r  coh 

and  the  frequency  selective  bandwidth  are  independent  of  the  propagation 


The  decorrelation  distance  of  the  signal  as  it  emerges  from  the 
delta  layer  is  approximately  equal  to  Under  strong  scattering  con¬ 

ditions  where  »  1,  this  distance  will  be  much  smaller  than  the  elec¬ 
tron  density  fluctuation  scale  size  and  the  quadratic  phase  structure 
approximation  conditions  are  met.  Conversely,  under  weak  scattering  con¬ 
ditions  where  -  1  the  quadratic  phase  approximation  will  give  inaccu¬ 
rate  results  for  the  mutual  coherence  function. 


Consider  now  the  mutual  coherence  function  for  two-positions, 
one-frequency  and  one-time: 

Mx.y.O.O)  *  exp  [-(x/zq)2  -  ( 5y/AQ) 2 ]  .  (2-80) 

The  distance  sl0  Is  then  seen  to  be  the  1/e  point  of  r^x, 0,0,0)  and  is 
therefore  the  x  direction  decorrelation  distance.  The  distance  i0/ 5  Is 
the  1/e  point  of  r^O.y.O.O)  and  is  therefore  the  y  direction  decorrela¬ 
tion  distance.  Because  6  is  less  than  or  equal  to  unity,  the  y  direction 
decorrelation  distance  is  always  greater  than  or  equal  to  l0. 

2.8  GENERALIZED  POWER  SPECTRAL  DENSITY. 


The  generalized  power  spectral  density  S(i^,T,ti>0)  of  the  signal 
incident  on  the  plane  of  the  receiver  Is  the  Fourier  transform  of  the 
mutual  coherence  function: 


$(fc,,T,u>0)  *  (2tr ) - 1  //  d2P  /  du>d  /  dt(jr1(p,a><j,td)  x 

’  -00  —OS  —OS 


exp  [-i(^*p-u)dT+a>0td)] 


(2-81) 


where  delay  t  is  the  Fourier  transform  of  frequency  difference  and 
Doppler  frequency  u>n  is  the  Fourier  transform  of  time  difference  t  ..  The 
quantity  S( i^,x ,ujp)d2K^dTd(Dp  is  equal  to  the  mean  signal  power  arriving 


within  the  angular  Interval  ^  to  +  d2i^;  with  delays  relative  to  a 
nominal  propagation  time  in  the  interval  t  to  t  +  dr;  and  with  Doppler 
frequencies  in  the  Interval  u>p  to  Up  +  dwp.  The  delay  spectrum  of  the 
GPSO  Is  a  consequence  of  the  fact  that  different  frequencies  propagate 
differently  through  the  ionosphere  with  some  frequencies  arriving  early 
and  some  frequencies  arriving  late.  The  Importance  of  this  effect  depends 
on  the  ratio  of  the  signal  bandwidth  to  the  frequency  selective  bandwidth. 


The  Indicated  integrals  in  Equation  2-81  can  all  be  done  in 
closed  form.  If  the  Kx  integral  is  done  first,  the  distance  offset 
results  In  a  term  with  the  form  exp{  iK^t^)  times  a  function  of  Kx .  The 
Fourier  transform  from  t^  to  Up  results  In  the  delta  function  2ir6p(up- 
Kxvg).  The  GPSD  is  then  the  product  of  a  Doppler  term  and  an  angular- 
delay  term  of  the  form 

S(Kx,Ky,T,Wp)  *  Sp(a1p)S(Kx,Ky,T)  .  (2-82) 

The  angular-delay  part  of  the  GPSD  Is 

S(Kx,Ky,T)  -  ;-/?) 1/2  <■»„„#-'  exp  t-(K**«-*K*)*|/*I  X 

exp  (-(a!/2)[»cot]t-A(K^K2)t2/4]2|  (2-83) 


where  the  delay  parameter  a  is  defined  to  be 


V0 
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(2-84) 


When  It  Is  recalled  that  Is  the  Rayleigh  phase  variance,  the  param¬ 
eter  a  Is  equal  to  1/C x  where  the  notation  Cj  Is  used  by  Wlttwer  (1980, 
1982,  1986).  The  components  of  are  related  to  the  scattering  angles 

e  and  8  about  the  x  and  y  axis  respectively  by  the  relations 
x  y 


Kx  =  2*  s1n(ex)/x 


K  ■  2«  s1n(0  )/X 

y  y 


(2-85) 
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For  the  frozen-in  approximation  and  assuming  that  the  effective 
velocity  is  along  the  x  direction,  the  Doppler  spectrum  is 


SD<“0>  =  !VfI"oV(iV 


(2-86) 


The  quantity  r0  is  the  decorrelation  time  of  the  signal  at  the 
receiver.  It  Is  formally  defined  here  to  be  the  time  required  for  the 
random  diffraction  pattern  to  drift  one  decorrelation  distance: 

lo  *  Vo  •  <2'87> 

However,  as  is  the  case  for  the  decorrelation  distance  and  the  coherence 
bandwidth,  the  decorrelation  time  is  calculated  using  the  more  general 
formalism  of  Wittwer. 


The  significance  of  the  parameters  a,  Zq  and  u>cQh  that  appear  in 
the  angular-delay  part  of  th  GPSD  will  be  discussed  in  the  next 
subsections. 

2.8.1  Delay  Spread  and  o. 

The  delay  spread  S(t)  of  the  signal  energy  arriving  at  a  fixed 
anqle,  given  by  the  second  exponential  term  in  Equation  2-83,  has  the 
Gaussian  form 

S( x )  =  («-coh/^)  exp  ( "a2wc0h( T"t p) 2/2 1  (2-88) 

where  tp  is  the  additional  propagation  time  for  signals  arriving  at  the 

angles  K  and  K  .  To  see  this,  the  expression  for  t  is  expanded  usinq 
x  y  p 

the  definitions  for  l  and  u>cQh  giving 

tp  *  (9x+ey)(zt+zr)(zr/-t)/2c  .  (2-89) 

The  geometry  of  the  scattering  in  one  plane  is  shown  in  Fioure  2-3.  It  is 
clear  that  for  small  angle  scattering,  the  anale  ©t  is  related  to  the 
receiver  scattering  anqle  by 
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Figure  2-3.  Scattering  geometry 


6t  =  (2r/zt)6r  *  (2_9n) 

The  difference  d  between  the  1 ine-of -sight  distance  and  the  scattered  path 
length,  for  small  angle  scattering,  is  given  by 

d  =  (zJ:+zV)1/2  +  (zt+zt9t) 1/2  “  ^Zr+Zt)  "  (zr/zt)(zr+2t)9r/2  (2_91) 

for  scattering  only  in  one  plane.  When  scattering  about  both  the  x  and  y 

axes  is  taken  into  account,  the  total  path  difference  is  given  by  Equation 

2-90  with  e2  replaced  by  02  +  Q2.  The  additional  time  required  for  the 
r  x  y 

signal  to  propagate  along  the  scattered  path  is  just  d/c  which  is  equal  to 


For  a  given  value  of  the  coherence  bandwidth,  larger  values  of  a 
mean  that  the  signal  energy  arrives  with  a  narrower  distribution  in  delay 
about  the  time  tp.  The  delay  parameter  a  is  then  a  measure  of  the  rela¬ 
tive  importance  of  diffraction  or  scattering  and  refraction  with  larqe 
.values  of  a  indicating  strong  scattering  effects  and  small  values  indi¬ 
cating  weak  scattering  or  refractive  effects  (Knepp  1982).  The  stronq 
scattering  limit  then  requires  that  the  values  of  2  be  large. 


2.8.2 


Frequency  Selective  Bandwidth  and  « 
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The  frequency  selective  bandwidth  f0  is  an  important  measure 
of  the  effects  of  scintillation  on  the  propagation  of  wide  bandwidth 
signals.  This  quantity  is  defined  as 

f0  -  1/(2ttot)  (2-92) 

where  orT,  the  time  delay  jitter  of  the  received  signal,  is 

02  =  <t2>  _  <t>2  _  (2-93) 

These  delay  moments  can  be  calculated  directly  from  the  angular-delay  part 
of  the  GPSD  using  the  definition 

OD  00 

P0<Tn>  =  (2*)"2  //  d2i<j^  /  dT  Tn  S(Kpt)  .  (2-94) 

-00  -  00 

It  is  easy  to  show  that  the  mean  received  power  P0  is  equal  to  unity. 

The  first  and  second  moments  are  also  obtainable  in  closed  form  giving  the 
relationship  between  w  h  ano  the  frequency  selective  bandwidth: 

u  u  ■  2*f  (1+l/a2)1/2  .  (2-95) 

coh  o 

This  expression  is  valid  only  under  the  quadratic  phase  structure  approxi¬ 
mation  because  Equation  2-94  has  been  evaluated  using  the  GPSD  calculated 
with  this  approximation.  Veh  and  Liu  (1977)  have  calculated  an  expression 
for  the  time  delay  jitter  keeping  both  the  second  order  and  the  fourth 
order  terms  In  the  expansion  of  A(p).  This  results  in  having  more  terms 
in  the  expression  for  the  time  delay  jitter.  However,  these  additional 
terms  will  be  significant  only  when  the  quadratic  approximation  for  A(p) 
is  invalid  and  therefore  only  when  the  expression  for  the  GPSD  is  also 
inval id . 
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The  1  +  1/a2  term  in  the  expression  for  the  coherence  bandwidth 
represents  the  relative  contributions  to  the  time  delay  jitter  from 
diffraction  (1)  and  refraction  (1/a2).  In  the  limit  that  a  is  large,  the 
time  delay  jitter  is  determined  by  diffractive  effects  alone  which  should 
be  the  case  under  strong  scattering  conditions. 

2.8.3  Angle-of -Arrival  Fluctuations  and  l  • 

o 

A  key  parameter  in  determining  the  severity  of  antenna  filtering 
effects  is  the  standard  deviation  o9  of  the  anqle-of-arrival  fluctua¬ 
tions  of  the  electric  field  incident  on  the  antenna.  It  is  clear  that  for 
anisotropic  scattering,  the  values  of  cr9  for  scattering  about  the  x  and 
y  axes  will  differ.  The  variance  of  the  angle-of-arrival  fluctuations 
about  the  5  direction  is  defined  as 

00  uo 

=  (2")'2  If  d2Kj_  /  dT  (XKc/2ir)2S(i(pT)  (2-96) 

•  00  .00 

under  the  small  angle  scattering  approximation  that  is  required  for  the 
6PSD  to  be  valid.  The  standard  deviations  of  the  angle-of-arri val 
fluctuations  about  the  x  and  y  axes  are 

o0x  -  */(✓!  ir£0)  (2-97) 

and 


°9y  s  x/(/2  *y5)  •  ( 2-98) 

For  the  small  angle  scattering  to  be  valid,  jQx  (which  is  the 
larger  of  the  two)  must  be  small  relative  to  1  radian.  Thus  the 
decorrelation  distance  i0  must  be  approximately  equal  to  or  greater  than 
wave  length  A.  The  small  angle  approximation  has  been  used 
throughout  the  derivation  of  the  GPSD,  starting  with  the  parabolic  wave 


equation.  The  resulting  expression  for  the  GPSD,  however,  does  not 
exhibit  singular  behavior  when  the  angle-of-arrival  fluctuations  become 
large  and  thus  the  small  angle  approximation  is  quite  robust. 

2.8.4  An  Isotropic  Example. 

When  the  penetration  angle  is  zero,  6  equals  unit  and  the 
scattering  is  isotropic  about  the  1 ine-of-sight.  The  one-dimensional 
generalized  power  spectral  density  is  then  given  by  the  integral 

09 

S(K,x)  =  (2*)-1  /  S(K,K  ,t)  dK 

*00  J 

■  («,'Sohy*,/v/1>  e*p  E(i/2“2)-<“cohT]  * 

Ff[l*o2(K2‘o/4-“c0h,)]/(2‘/!“)]t  (2-99) 

where  the  function  F  is  defined  as 

oo 

F(z)  *  exp  (-z2)  /  exp  (-t4-2t2z)  dt  .  (2-100) 

*  00 

The  F  function  has  been  approximated  by  Wittwer  (1980)  usinq  a  polynomial 
expansion  that  is  accurate  to  within  one  percent.  This  function  may  also 
be  written  in  terms  of  Ki/4  and  Ij/t*  Bessel  functions  (Knepp  1982). 

A  three  dimensional  plot  of  the  isotropic  one-dimension  GPSD 
with  a  equal  to  4  is  shown  in  Figure  2-4.  This  plot  shows  the  mean 
received  power  as  a  function  of  angle  K  and  delay  t.  The  vertical  axis  is 
linear  with  arbitrary  units. 

It  can  be  seen  that  the  power  arriving  at  large  angles  is  also 
the  power  arriving  at  long  delays.  The  power  arriving  at  long  delays  thus 
has  higher  spatial  frequency  components  than  power  arriving  at  short 
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Figure  2-4.  Generalized  power  spectral . density. 

delays.  Under  the  frozen-in  approximation  where  the  ionosphere  is  modeled 
as  a  rigid  structure  drifting  across  the  1  ine-of-sight,  these  higher  spa¬ 
tial  frequency  components  correspond  to  higher  Doppler  frequency  compo¬ 
nents.  The  signal  energy  arriving  at  long  delays  then  varies  more  rapidly 
in  time  than  the  signal  energy  arriving  at  short  delays  under  the  frozen- 
in  approximation. 

2.9  TURBULENT  APPROXIMATION. 

Under  conditions  before  striations  have  formed  in  the  ionosphere 
or  when  there  are  multiple  scattering  layers  with  different  velocities  in' 
the  line-of-sight,  the  frozen-in  approximation  may  provide  a  poor  model  of 
the  temporal  variations  of  the  received  signal.  An  alternative  to  the 
frozen-in  approximation  is  the  turbulent  approximation  where  the  temporal 


variations  and  the  spatial  variations  of  the  received  signal  are  independ¬ 
ent.  The  mathematical  form  of  this  approximation  is  that  the  two-posi¬ 
tion,  two-frequency,  two-time  mutual  coherence  function  is  separable  into 
a  product  of  a  spatial -frequency  term  and  a  time  term: 

rx(p»"(J.td)  «  r(p,u>d)r(td)  .  (2-101) 

The  coherence  function  r(p,u>d)  is  given  by  Equation  2-75  with  the  time 
difference  tj  set  to  zero.  After  performing  the  necessary  Fourier 
transforms,  the  6PSD  for  this  model  has  the  form 

S(i<pT,a)0)  =  Sq(  uq)  S(Kj^,  t)  (2-102) 

where  for  this  model  the  Doppler  spectrum  Sq(wq)  is  a  function  only  of  the 
Doppler  frequency.  Thus  Sp(a)p)  does  not  couple  angles  and  Doppler  fre¬ 
quencies  or,  equivalently,  positions  and  times  as  is  the  case  with  the 
frozen-in  approximation.  The  angle-delay  term  in  Equation  2-102  is  then 
qiven  by  Equation  2-S3. 

Under  the  frozen- in  approximation,  the  temporal  variations  of 
the  signal  at  long  delay  are  more  rapid  than  the  temporal  variations  of 
the  signal  arriving  at  short  delays.  Under  the  turbulent  approximation, 
the  Doppler  spectrum  is  independent  of  delay  so  the  temporal  variations  at 
all  delays  have  the  same  rate. 


In  order  to  specify  the  Doppler  spectrum  under  the  turbulent 
approximation,  the  two  position,  two  time  mutual  coherence  function  of  the 
electron  density  fluctuations  is  needed.  Detailed  descriptions  of  the 
spatial  variations  of  the  electron  density  fluctuations  are  currently 
available  (Wittwer  1986).  However,  detailed  information  on  the  temporal 
variations  of  the  electron  density  fluctuations  in  the  ionosphere  is  not 
currently  available.  An  f*4  Doppler  spectrum  will  therefore  be  assumed 
for  implementation  convenience  in  Section  5  where  channel  simulation 
techniques  are  discussed. 
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2.10 


IMPULSE  RESPONSE  FUNCTION  AND  ANTENNA  EFFECTS. 


The  channel  impulse  response  function  of  the  signal  Incident  on 
the  plane  of  the  receiver  and  the  impulse  response  function  of  the  signal 
at  the  output  of  an  aperture  antenna  will  be  discussed  in  this  subsection. 

2.10.1  Channel  Impulse  Response  Function. 

Consider  a  solution  U( p,zr,w,t)  to  the  parabolic  wave  equation 
in  the  plane  of  the  receiver.  This  represents  the  random  effects  due  to 
the  fluctuating  ionosphere  on  the  incident  electric  field  at  position  p 
and  time  t  from  a  transmitted  monochromatic  wave  with  angular  frequency 
w.  The  channel  impulse  response  function  of  the  signal  in  the  receiver 
plane  is  (Knepp  and  Wittwer  1984) 

00 

h(p,T,t)  =  (2tt) - 1  /  U(p,2r,«+«0,t)  exp  [ie(w)+iu)T]  dw  (2-103) 

-CO 

where  u>Q  is  the  carrier  anaular  frequency  and  6(a)  is  the  dispersive  con¬ 
tribution  to  the  impulse  response  function  due  to  the  mean  ionization. 

The  term  exp  [  i e ( u>)  ]  is  the  transfer  function  of  a  smooth  ionized  plasma 
and  is  equal  to  the  exponential  term  in  Equation  2-16.  Thus  6(w)  is 

Zr__  2r 

e(w)  =  -  /  k( z * ) dz '  =  -  (»/c)  /  [l-'/(z'  )/<*)2] 1/2  dz '  .  (2-104) 

-*t  -*t  P 

Because  the  smooth  plasma  or  dispersion  effects  represented  by  exp[i6(u>)] 

and  the  fluctuating  plasma  effects  represented  by  U(p,zr,w,t)  appear  as 
the  product,  it  is  convenient  to  separate  these  effects.  The  dispersive 
effects  will  be  considered  in  Section  2.10.2. 
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If  the  transmitted  signal  is  a  modulated  waveform  m(t)  then  the 
signal  complex  voltage  incident  on  the  plane  of  the  receiver  is  the  convo¬ 
lution  of  the  transmitted  modulation  and  the  channel  impulse  response 
function: 


e(p,t)  =  /  m(t-T-t  )h(p,x,t)  dt  (2-105) 

_oo  ” 

where  tp  is  the  nominal  propagation  time.  If  the  delay  spread  of  the 
impulse  response  function  is  larger  than  the  symbol  period  of  m(t),  then 
the  convolution  will  encompass  multiple  symbols  with  intersymbol  inter¬ 
ference  as  a  result.  It  is  also  clear  from  this  equation  that  signal 
energy  arriving  at  longer  delays  corresponds  to  symbols  transmitted  at 
earlier  times. 

2.10.2  Dispersive  Effects. 

When  the  dispersive  term  F(u>)  is  expanded  in  a  Taylor  series 
about  the  carrier  radian  frequency,  the  result  is 

9  ( (u)  =  0(w  )  -  (u>-u>  )  9  *  ( u>  )  +  (w-u>  )20"(u)  )/2  +  ...  (2-106) 

where  the  first  three  coefficients  in  the  expansion  are 

zr  _7 

9(o>  )  =  -((o  /c)  /  [l-o)  (z')/u>2]1/2  dz'  (2-107a) 

zr  _2 

6'(u)  )  =  (1/c)  j  [  z '  )/o)2]"1/2  dz'  ( 2-  107b ) 

0  -z>  p 

v 

Zr 

If" ( o>o )  =  -( 1/c )  /  (^p/oj^)[l-o >p(z' )/o)2]*1/2  dz'  .  (2-107c) 
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These  equations  may  be  expanded  using  the  assumption  that  the  carrier  fre 
quency  is  much  larger  than  the  plasma  frequency.  The  first  three  coeffi¬ 
cients  then  reduce  to 


8(ioq)  •  -2itR/A  + 

(2-108a) 

?'(«!))■  R/c  +  \2r  N_/(2irc) 
o  e  T 

(2-108b) 

?"(«  )  -  -A3reNT/(2it2C2) 

(2-108C) 

where  the  free  space  range  R  and  total  electron  content  (TEC) 

Ny  are 

and 

R  =  z.  +  z 
t  r 

(2-109) 

HJ  =  V  <"e(z')>  dz*  . 

"zt 

(2-110) 

The  first  terms  in  Equations  2-108a  and  2-108b  are  simply  the  free  space 
phase  shift  and  propagation  time  which  are  proportional  to  the  line-of- 
sight  distance  R.  The  terms  proportional  to  Nj  in  Equations  2-108 
represent  the  phase  shift,  group  delay  and  dispersion  due  to  the  mean 

ionization. 

The 

TEC  dynamics 

Doppler  shifc  fp  of  the  incident  signal  due  to 

is 

range  and 

fQ  -  (2it )_1  de(w  )/dt  -  -  2nR/ A  +  ArgNT  . 

(2-111) 

Note  that  increasing  TEC  (positive  Nj)  increases  both  the  propagation 
time  and  ‘.he  Doppler  shift  whereas  increasing  R  (positive  R)  increases  the 
propagation  time  but  decreases  the  Ooppler  shift. 


2.10.3  Antenna  Aperture  Infects. 


The  voltage  at  th>  a  of  an  aperture  antenna  is  the  spatial 
convolution  of  the  incident  voltage  and  the  aperture  weiqhting  function. 
The  received  voltage  for  an  antenna  located  at  pQ  ana  pointing  in  the  i<0 
direction  is  then  given  by  (Knepp  1983a) 

00 

uA(P0»“.t)  *  //  U(p\zr,w,t)Ar(P0-P')  exp  ( iKQ*P ' )  d2p'  (2-112) 

•  oo 

where  the  subscript  A  denotes  the  voltage  at  the  output  of  the  antenna. 

The  z  dependence  of  U/\  has  been  suppressed  because  it  is  understood  that 
chis  voltage  is  at  the  receiver  plane.  It  is  assumed  that  the  aperture 
weighting  function  of  the  receiver  Ar(p)  is  independent  of  frequency. 

This  is  generally  true  for  a  range  of  frequencies  about  the  carrier  fre- 
auency  that  is  larger  than  the  signal  bandwidth. 

In  order  to  relate  the  6PS0  of  U/\  to  the  fiPSD  of  the  incident 
signal,  the  two-position,  two-frequency,  two-time  mutual  coherence  func¬ 
tion  of  'J/\  is  reauired.  The  mutual  coherence  function  of  the  signal  out 
of  the  antenna  is 

00  00 

//  d2p ’  //  d2p"  <U(p',zr,u)1,t1)U*(p,,,zr,u)2,t2)>  * 

_,U0  .00 

Ar(pi-p')Ar*(p2-p")  exp  [  1I<0*  (p'-p")  ]  .  (2-113) 

For  statistically  stationary  processes,  the  expectation  of  UAUA*  must  be 
a  function  only  of  the  differences  -  oj  -  02,  =  oi x  -  w2,  and  t^  = 

t i - 1 2  and  the  expectation  of  UU*  in  the  integrand  must  be  a  function  only 
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of  the  differences  p'  -  p",  o>d,  and  t^.  The  aperture  weighting  function 
may  be  written  in  terms  of  the  angular  beam  profile  using  the  Fourier 
transform  relationship 

OD 

Ar(p)  =  (2ir)-2  //  Ap(^)  exp  (1^*p)  d2^  (2-114) 

Upon  substituting  this  equation  for  both  aperture  weighting  functions  In 
the  expression  for  the  mutual  coherence  function,  changing  variables  from 
p'  to  p  a  p‘  -  p",  and  changing  the  order  of  integration,  Equation  2-113 
becomes 

00 

VVVV  '  (2ir)'4  H  r^*zr’“d’td)  exp  <V>  x 

.00 

00  « 

//  d2K  Ar*(i<)  exp  [iK*(Pi-p)]  //  d2K*  Ar(K')  exp  [-1K ’ -p2]  * 

.00  .« 


//  d2p"  exp  [ip" • (K1 -K) ]  .  (2-115) 

-W 

The  last  Integral  In  this  expression  is  equal  to  (2")2Sp(K'-K)  *nd  the  K’ 
integral  may  be  performed  directly.  Another  change  in  the  order  of 
integration  results  In 

oo 

VVVV  s  M  Gr(^  exp  (1^'^d)  * 

•  «o 

00 

//  d2P  r(p.zr.«d.td)  exp  (-1p*(£-i<0)]  .  (2-116) 

.00 

The  quantity 

Gr(K)  -  Ar*(K)Ar(K)  (2-117) 


is  the  power  beam  profile  of  the  receiving  antenna. 


The  mutual  coherence  function  r(p,zr,u)d,fd)  of  the  signal  Inci¬ 
dent  on  the  plane  of  the  antenna  that  appears  In  the  second  integral  of 
Equation  2-116  Is  the  product  of  the  free  space  term  r0  (Equation  2-43) 
and  the  stochastic  term  r1#  The  free  space  term  may  be  pulled  out  of  the 
second  Integral  if  It  is  assumed  not  to  vary  over  the  face  of  the  antenna. 
This  is  equivalent  to  assuming  that  any  deviations  from  a  plane  wave  In 
the  Incident  signal  are  due  to  scattering  effects  In  the  Ionosphere  and 
are  not  due  to  geometrical  effects.  After  the  free  space  term  is  pulled 
out  of  the  integral,  rA  may  be  assumed  to  represent  only  the  stochastic 
fluctuations  of  the  received  signal. 

Now  the  GPSD  of  the  signal  out  of  the  antenna  may  be  computed  by 
taking  the  appropriate  Fourier  transforms  (see  Equation  2-81)  from  Pd,  u>d, 
and  td  to  K^,  t,  and  u>p  respectively.  After  performing  the  u»d  to  t  trans¬ 
form,  the  angular-delay  part  of  the  GPSD  at  the  antenna  output  will  be 

oo 

sA(KpT)  s  (2")“2  II  d2i<  Gr(K)  x 

00  00 

//  d2p  r1(p,zr,T)  exp  [-i(K-i<0)*p]  II  d2Pd  exp  [i(K-K,  )*Pd]  .  (2-118) 

—00  M  00 

The  last  Integral  In  this  equation  Is  just  ( 2tt ) 2 5p( i<-^j_)  and  the  middle 
Inteciral  gives  the  angular-delay  GPSD  of  the  incident  voltaoe  S(K-£0,t). 
The  GPSD  of  the  signal  out  of  the  antenna  is  then 

SA(KpT)  -  Gr(^)S(Kj_40,T)  .  (2-119) 

The  effect  of  an  antenna,  as  should  be  expected,  Is  to  modify  with  the 
beam  profile  the  mean  incident  power  as  a  function  of  angle.  This  result 
will  be  used  throughout  the  rest  of  this  report. 
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SECTION  3 

ANTENNA  FILTERING  EFFECTS 


An  antenna  beam  profile  acts  as  an  angular  filter  of  the  re¬ 
ceived  signal  energy.  Because  of  this,  the  mean  power,  decorrelation  dis¬ 
tances,  and  frequency  selective  bandwidth  of  the  signal  at  the  output  of 
an  antenna  are  different  than  those  of  the  incident  signal.  The  reduction 
In  mean  power  Is  a  direct  consequence  of  the  attenuation,  due  to  the  beam 
profile,  of  energy  arriving  at  large  angles-of-arrlval  relative  to  the 
peak  of  the  beam.  The  energy  arriving  at  large  angles  is  also  the  energy 
arriving  at  long  delays.  The  frequency  selective  bandwidth  Is  an  inverse 
measure  of  the  delay  spread  of  the  signal  energy.  Hence  the  frequency 
selective  bandwidth  of  the  signal  out  of  an  antenna  Is  larger  than  the 
frequency  selective  bandwidth  of  the  signal  incident  on  the  antenna.  The 
decorrelation  distance  is  an  Inverse  measure  of  the  angle-of-arri val  fluc¬ 
tuations  of  the  signal  energy.  The  effect  of  an  antenna  Is  to  reduce  the 
angular  spread  of  the  signal  and  thus  to  Increase  the  decorrelation  dis¬ 
tances  of  the  signal  out  of  the  antenna  relative  to  the  decorrelation 
distances  of  the  Incident  signal. 

The  effects  of  aperture  antennas  with  arbitrary  beamwidths  will 
be  considered  in  this  section.  The  antenna  beam  profiles  for  uniformly 
weighted  circular  or  rectangular  apertures  and  for  Gaussian  apertures  are 
described  in  Section  3.1.  The  filtering  equations  for  mean  power,  spatial 
and  temporal  decorrelation  properties,  and  frequency  selective  bandwidth 
of  the  signal  out  of  a  Gaussian  antenna  with  arbitrary  beamwidths  are 
derived  In  Section  3.2.  Then  a  comparison  of  the  antenna  filtering 
effects  of  a  Gaussian  antenna  and  a  uniformly  weighted  circular  antenna  is 
made  in  Section  3.3  to  assess  the  accuracy  of  the  Gaussian  approximation. 
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3.1 


ANTENNA  DESCRIPTIONS. 


The  coordinate  systems  for  the  propagation  and  the  antenna  are 
shown  in  Figure  3-1.  The  z  axis  is  along  the  line-of-slght  between  the 
transmitter  and  the  receiver  and  the  direction  with  the  minimum  decorrela¬ 
tion  distance  Is  along  the  x  axis.  The  antenna  u-v  coordinate  system.  In 
this  case  for  a  rectangular  antenna,  is  In  the  x-y  plane.  The  rotation 
angle  x  Is  the  angle  between  the  x  axis  and  the  u  axis. 

In  general,  the  face  of  the  antenna  will  not  lie  in  the  x-y 
plane.  However,  for  satellite  communication  links  It  is  usual  for  the 
satellite  position  to  be  known  accurately  through  ephermeris  data  and  for 
an  antenna  to  be  dedicated  to  a  single  link.  Thus  the  rotation  angles  of 
the  antenna  out  of  the  x-y  plane  will  be  small  in  most  cases  and  the 
cosine  squared  effects  of  these  angles  may  be  ignored.  If  these  rotation 


Figure  3-1.  Propagation  and  antenna  coordinate  systems. 
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angles  are  not  small,  then  in  the  developments  that  follow  the  antenna 
size  is  the  projection  in  the  x-y  plane  and  the  antenna  beamwldths  must 
include  the  effects  of  beam  broadening  as  the  beam  is  pointed  away  from 
boresight. 


3.1.1  6auss1an  Bean  Profiles. 

The  antenna  beam  profile  is  assumed  to  be  separable  in  the  u-v 
coordinate  system.  The  Gaussian  antenna  beam  profile  is  then 

G(KU,KV)  ■  exp  (-»*K*-°^)  .  (3-1) 

The  peak  gain  G(0,0)  has  been  set  to  unity  because  this  value  is  usually 
included  in  the  calculation  of  the  mean  received  power.  For  either  the  u 
or  the  v  direction,  the  antenna  beam  profile  can  also  be  written  as 

G ( 0  )  -  exp  [-(In2)(20/0O)2]  (3-2) 

where  0O  is  the  3  dB  beamwidth  (i.e.  full  width  at  half  maximum). 

Equating  these  two  profiles  gives 

a2  =  (ln2)X2/(n2025)  (3-3) 

where  0^  is  the  3  dB  beamwidth  in  the  £  =  u  or  the  £  a  v  direction. 

The  beam  profile  is  required  in  the  x-y  coordinate  system  of  the 
propagation.  Using  the  transformation 

K.  «  Kcosx  +  K  sinx 

y  (3-4) 

K  *  -  K  sinx  +  K  cosx 

v  a  y 

gives 
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G(Kx,Ky)  =  exp  [-(a*cos2x+aJsin2x)K2]  * 

exp  [-2(a2-aJ)sinxcosxKxKy] 

exp  [-(a2sin2x+a2cos2x)K2]  .  (3-5) 

The  coefficients  and  will  be  related  to  the  sizes  of  uniformly 
weighted  apertures  in  the  next  subsections.  If  a  and  ay  are  equal,  then 
the  Gaussian  beam  profile  is  independent  of  the  rotation  angle  and  is 
referred  to  as  isotropic  in  subsequent  sections. 


3.1.2  Uniformly  Weighted  Circular  Apertures. 


For  a  uniformly  weighted  circular  aperture,  the  aperture 
weighting  function  without  pointing  errors  is 


V5> 


1/a  if  | p j  <  D/2, 
.  0  otherwise 


(3-6) 


where  D  is  the  diameter  of  the  circular  aperture  and  where  x  =  r  for  a 
receiving  antenna  and  x  =  t  for  a  transmitting  antenna.  The  value  of  a  is 
chosen  so  that  the  peak  antenna  gain  G(0)  is  unity.  The  voltage  antenna 
gain  pattern  Ax(i<)  is  related  to  the  aperture  distribution  function  by  the 
Fourier  transform 


Ax(i<)  *  //  exp  (-ii<*p)Ax(p)d2p  .  (3-7) 

«oo 

The  coordinate  system  for  this  transformation  is  shown  in  Figure  3-2.  The 
antenna  is  pointed  along  the  1 i ne- of -sight  in  the  z  direction  and  p  is 
chosen  to  be  along  the  x  axis.  The  beam  profile  for  a  circular  antenna 
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Figure  3-2.  Circular  antenna  coordinate  system. 

will  then  be  a  function  only  of  the  elevation  angle  0.  In  this  coordinate 
system,  the  dot  product  is 


£«p  a  kr  sine  s  1  n<t> 


(3-8) 


where  k  3  |£|  *  2it/x  and  r  =  |pj.  The  Fourier  transform  then  becomes 

D/2  » 

A  (0)  =  a'2  /  dr  r  /  d$  exp(-ikr  sine  sin<j>)  .  (3-9) 

x  o  -it 

Performing  the  indicated  integrals  results  in  the  well  known  form  for  the 
power  beam  prof 1 le 


G  (0)  3  A  ( 0 ) A  *(0)  =  4JL2i1Tii/.x)si.n9I 
x  4  x  U(D/x)sine]2 


(3-10) 


f'j  :  -■  v  f  /vv  v  ./■  v 


where  Ji  is  the  Bessel  function  and  where  6(0)  =  1  when  the  value  of  a  is 
chosen  to  be  the  area  of  the  aperture  ttD2/4.  The  factor  of  4  in  the 
expression  for  Gx(8)  is  required  because  the  limit  of  Ji2U)/52  as  % 
approaches  zero  is  1/2.  The  3  dB  full  width  at  half  maximum  beamwidth 
[i.e.  G(0o/2)  =  1/2]  is  given  by  solving  the  equation 

4Ji2(5)/C2  *  1/2  (3-11) 

with  the  result  5  =  1.616340.  Assuming  that  the  beamwidth  is  small  so 

sin(e  /2)  ■  0  /2,  the  beamwidth  in  terms  of  the  diameter  D  is 
o  o 

0Q  =  1.02899  (X/D)  radians  .  (3-12) 

If  this  beam  profile  is  approximated  by  a  Gaussian  profile  with  the  same  3 
dB  beamwidth,  then  the  a  coefficients  that  appear  in  Equation  3-1  are 

°12  =  =  ( ln2)D2/(l. 02899 it)2  .  (3-13) 

3.1.3  Uniformly  Weighted  Rectangular  Apertures. 


For  uniformly  weighted  rectangular  apertures,  the  aperture 
weighting  function  without  pointing  errors  is 

|  1/a  if  | u |  <  Du/2  and  |v|  <  Dy/2 
Ax(u,v)  =  (3-14) 

I  0  otherwise 

where  Du  and  Dv  are  the  lengths  of  the  aperture  in  the  u  and  v  directions 
respectively.  In  this  case,  the  Fourier  transform  indicated  in  Equation 
3-7  can  be  done  directly  with  the  result 


G  (8  ,0  )  -  G  (0  )G  (0  ) 
x'u’v'  x '  u  x  v ' 


(3-15) 


where  G  (6,)  for  5  =  u  or  5  =  v  has  the  familiar  sin2(x)/x2  form 

X  s 


sin2[Tr(0{./X)sin3{.  ] 
G  (0.)  = _ i _ 

U(D./x)sin0£  ]2 


(3-16) 
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The  normalization  a  of  the  weighting,  function  is  just  the  area  D^D  of  the 
rectangular  antenna.  The  3  dB  beamwidth  is  given  by  solving  the  equation 

sin2U)/52  *  1/2  (3-17) 

with  the  result  5  =  1.391557.  The  e  direction  beamwidth  in  terms  of  the 
antenna  size  0^  is  then 

0Q5  =  0.885893  (x/0  )  radians  (3-18) 

assuming  that  sin(0o^/2)  -  0^/2.  If  this  beam  profile  is  approximated 
by  a  Gaussian  profile  with  the  same  3  dB  beamwidth,  the  a  coefficients 
that  appear  in  Equation  3-1  are 

a2  =  (ln2)D2/(0.885893*)2  .  (3-19) 

S* 

3.2  FILTERING  EQUATIONS. 

The  filtering  equations  relate  the  statistics  of  the  signal  at 
the  outputs  of  one  or  more  antennas  to  the  statistics  of  the  signal  inci¬ 
dent  on  the  antennas.  The  statistics  that  are  considered  in  this  section 
are  the  mean  power,  frequency  selective  bandwidth,  decorrelation  dis¬ 
tances,  and  decorrelation  times  of  the  signal  out  of  an  antenna  and  the 
cross  correlation  of  the  signals  out  of  two  separate  antennas.  A  Gaussian 
beam  profile  is  used  because  this  approximation  leads  to  closed  form 
expressions  for  the  filtering  equations.  The  accuracy  of  this  approxima¬ 
tion  is  investigated  in  Section  3.3. 

The  power  impulse  response  function  G/\(x)  of  the  signal  out  of 
an  antenna  is  also  calculated  in  this  section  in  the  limit  that  the  param¬ 
eter  a  is  equal  to  infinity.  The  quantity  G^-Odt  is  equal  to  the  mean 
received  power  in  the  delay  interval  i  to  i  +  dt.  This  function  will  be 
used  in  Section  5  which  describes  channel  simulation  techniques. 
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The  severity  of  the  filterinq  effects  is  determined  by  the  rela¬ 
tive  size  of  the  standard  deviation  of  the  anqle-of-arrival  fluctuations 
aQ  and  the  antenna  beamwidths.  When  is  small  compared  to  the  antenna 
beamwidths,  the  signal  arrives  essentially  at  the  peak  of  the  beam  pro¬ 
file,  if  the  pointing  error  is  small,  and  the  filtering  effects  are  small. 
However,  if  o6  is  large  compared  to  the  beamwidths,  much  of  the  signal 
energy  arrives  at  angles  outside  of  the  main  lobe  of  the  antenna  beam  pro¬ 
file  and  the  filtering  effects  are  large.  Equivalently,  large  values  of 
the  ratio  ae/e0  correspond  to  situations  where  the  decorrelation  distance 
of  the  incident  signal  is  small  compared  to  the  antenna  size  and  the  inci¬ 
dent  electric  field  as  seen  by  the  aperture  is  no  longer  a  plane  wave.  In 
this  situation,  the  induced  voltages  across  the  face  of  the  aperture  do 
not  add  coherently  when  summed  together  by  the  antenna  with  a  loss  in 
signal  power  as  a  result. 

3.2.1  Orthogonal  1  zed  GPSD. 

The  delay  moments  of  the  signal  out  of  an  antenna  pointing  at 
the  transmitter  are  given  by 

oo  OO  00 

PA<Tn>  =  (2*)-3  //  d2^  Gr(KjJ  /  dT  t n  /  du>D  S(i^tTfw0)  (3-20) 

where  P/\  is  the  mean  power  of  the  signal  out  of  the  antenna  and  is  equal 
to  the  n  =  0  delay  moment.  The  Doppler  integral  is  just 

oo 

S(Kj_,t)  =  /  S{K^,T,u>D)du>D  (3-21) 

where  S(K^,t)  is  the  angular-delay  part  of  the  GPSD  of  the  incident  signal 
given  by  Equation  2-83. 


When  Equations  2-83,  2-119,  and  3-5  are  combined,  the  resulting 
GPSD  of  the  signal  out  of  the  antenna  in  the  x-y  coordinate  system  is 

SA(ltx*VT)  =  ('/2,l/2  ““cot/o5'1  exp  (-(a2/2)[“cohT'A(KxtKJ)£o/4!2)  * 

exP  (-42K2x«2/4]  exp  [-«Jk2iJ/452]  exp  (3-22) 

where  the  a  coefficients  are 

A2  =  1  +  4a2cos2x/t2  +  4a2sin2x/*2 
x  u  o  v  o 

A2  =  1  +  462a2Sin2x/t2  +  462a2C0S2x/^2 

Axy  =  S6(a2-a2)sinxcosx/t2  .  (3-23) 

The  cross  term  K  K  in  the  last  exponent  of  Equation  3-22  signi- 
x  y 

ficantly  complicates  the  evaluation  of  the  angular  integrals.  However, 

this  term  may  be  eliminated  by  performing  the  integration  in  a  rotated  p-q 

coordinate  system  where  the  GPSD  is  orthogonal  (i.e.  contains  only  K2  and 
,  * 

Ky  terms).  The  rotation  is  defined  by 

K  -•  K  cose  +  K  sine 

X  P  Q  (3-24) 

K  *  -K  sine  +  K  cose 

y  P  9 

where  dKpdKq  =  dKxdK^  and  Kp  +  Kq  =  K2  +  Ky'  ^he  an9^e  0*  ch°sen  so  that 

the  coefficient  of  is  zero,  is 

P  0 

e  =  (l/2)tan-1[25A2y/(52A2-A2)]  .  (3-25) 


In  the  rotated  coordinate  system,  the  GPSD  at  the  antenna  output  takes  the 
form 


3.2.2  Scattering  Loss. 

The  mean  power  of  the  signal  out  of  the  antenna  is  calculated 
from  Equation  3-28  with  n  =  0  and  is 

P.  =  1/A.A  =  l/l  (3-29) 

A  p  q  5 

where  L$  is  the  scattering  loss.  After  some  manipulations,  the  square 
of  the  scattering  loss  becomes 

l J  =  (l+462a2/t2)(l+4a2/t2)sin2x  +  ( l+4a2/ i2  )  ( l+4<52a2/t2  )cos2x  .  (3-30) 

S'  uo  vo  uo  vo 

The  terms  in  parenthesis  occur  frequently  in  the  filtering  equations. 

Thus  for  notational  convenience,  define 
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■  >  *  <81"2><V9ou> 


6uy  '  '  ♦  ««VX  =  1  *  <8'"2><  V9o/ 

6vx  *  1  +  4“^o  "  1  *  (81n2 ) ( c0>(/eov )2 

Gvy  *  1  ♦  4«V/*o  -  1  *  (81n2)(o9y/6ov)2  (3-31) 

where  equations  2-97  and  2-98  have  been  used  to  write  a  and  aq/5  in  terms 

of  the  variance  of  the  angle-of-arri val  fluctuations  o  and  o  and  where 

“x  ey 

Equation  3-3  has  been  used  to  write  and  in  terms  of  the  3  dB  beam- 

widths  e  and  0  .  With  these  definitions,  the  scattering  loss  takes  the 

cu  ov 

simpler  form 


=  [G  G  sin 
L  uy  vx 


Jx+G  G  cos2x]1/2 
A  ux  vy  J 


(3-32) 


The  scattering  loss  is  a  function  of  the  ratios  of  the  widths  of 
the  angie-of-arri val  fluctuations  oQ  and  the  antenna  beamwidths  0Q  with 
the  scattering  loss  approaching  unity  for  og  «  9q.  It  can  be  seen  that 
the  dependence  un  the  rotation  angle  disappears  if  the  scattering  is 
isotropic  (i.e.  6  =  I  and  oQx  =  o0y)  or  if  the  antenna  is  isotropic  (i.e. 

9ou  9ov^* 


Now  consider  the  case  of  an  isotropic  antenna  where  G^x  =  G^ 


When  the  scattering  is  also  isotropic  (penetration  angle 


and  G  =  G  . 
uy  vy 

=  O’),  G  =  G  and  the  scattering  loss  is  Lr  =  G  .  However  when  the 
’  ux  uy  S  ux 

penetration  angle  is  90*,  the  y  direction  decorrelation  distance  will  be  q 
( q- 1 5)  times  larger  than  the  x  direction  decorrelation  distance.  Then 


under  most  circumstances,  G  -  1  and  the  scattering  loss  is  L  »  /G  . 

uy  b  ux 

The  scattering  loss  is  therefore  a  sensitive  function  of  the  penetration 


angle. 


.  -  .V.V  , 


3.2.3 


Frequency  Selective  Bandwidth. 


The  frequency  selective  bandwidth  has  been  defined  in  terms  of 
the  time  delay  jitter  of  the  signal  in  Section  2.8.2.  The  required  delay 
moments  are 


u>  ,<t>  =  (A/2)  [a-2+62a-2] 
coh  '  '  L  p  q  J 

and 

“coh<T2>  =  [a'2+(3A2/4A4p)  ♦  (AV/2AV)  +  ( 3 A2 <5 **/ 4A4 )  ] 
The  time  delay  jitter  is  then  given  by 


T 


/2)(A-4+S4A“4) ] 


(3-33) 

(3-34) 

(3-35) 


The  a"2  term  in  the  expression  for  a2  repres<  ts  the  effects  of 

refraction  on  the  time  delay  jitter.  In  the  limit  that  diffraction 

effects  dominate  the  signal  fluctuations,  the  time  delay  spread  of  the 
signal  is  determined  by  the  spread  of  propagation  times  as  a  function  of 
the  angle-of-arri val  of  the  signal  with  the  largest  delays  corresponding 
to  the  largest  angles-of-arrl val .  The  effect  of  an  antenna  then  is  to 
preferentially  attenuate  the  signal  energy  at  large  angles  and  thereby  to 
reduce  the  delay  spread  of  the  signal  arid  to  Increase  the  frequency  selec¬ 
tive  bandwidth.  In  the  limit  of  a  large  angle-of-arrival  variance,  or 
equivalently  of  a  small  decorrelation  distance,  only  energy  arriving  near 
2ero  angle  is  received  at  the  output  of  the  antenna  and  the  time  delay 

jitter  of  the  output  signal  approaches  zero.  The  value  of  o  must  be  set 

large  enough  that  the  time  delay  j  er  of  the  signal  out  of  the  antenna 
is  not  limited  by  a  under  the  most  severely  disturbed  propagation  condi¬ 
tions.  The  expressions  below  for  thf-  frequency  selective  bandwidth  of  the 
signal  out  of  the  antenna  are  therefore  calculated  with  a  set  to  infinity. 


The  ratio  of  the  frequency  selective  bandwidth  of  the  signal  out 
of  the  antenna  to  that  of  the  Incident  signal  Is  just 

f./f  *  o  (without  antenna)/?  (with  antenna)  .  (3-36) 

A  0  T  T 

The  numerator  of  this  expression  may  be  obtained  from  Equation  3-35  with 
ap  and  set  to  unity.  The  resulting  equation  for  the  ratio  Is 

f./f  =  [(1+«-)/(1/a>V^)]1/2  .  (3-37) 

A  o  p  q 

After  some  manipulations,  this  expression  may  be  rewritten  as 

(1+64)1/2  L? 

f/f  » - £ -  .  (3-38) 

[(54G2  +G2  )cos2x  +  (G2  +64G2  )sin2x]1/2 
1  ux  vy  uy  vx  J 

3.2.4  Two-Position  Mutual  Coherence  Function. 

It  will  be  convenient  in  further  developments  to  have  the  two- 
position  mutual  coherence  function  of  the  signal  out  of  the  antenna  calcu¬ 
lated  in  the  x-y  coordinate  system.  The  coherence  function  is  defined  as 

OD  00  co 

rA(x,y)  =  (2*)-2  /  dKx  f  dK  (  dr  exp(1K  x+IK  v)Sa(K  K  t)  (3-39) 

rw  CO  —00  y  -00  y  y 

where  Equation  3-22  is  used  for  the  GPSO.  In  this  case,  the  integral  may 
be  done  in  closed  form  with  the  result 

rA(x,y)  «  PA  exp  {-[(Guysin2x+Gvycos2x)/L^]x2/z2)  * 
exp  {-((Guxcos2x+Gvxsin2x)/L2]62y2/t^}  * 

exp  {-[2(G|jx-Gvx)sinxcosx/L2]52xy/t2)  .  (3-40) 
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3.2.5 


Decorrelation  Distances  and  Time. 


The  x  direction  decorrelation  distance  of  the  signal  out  of  the 
antenna  Is  given  by  the  1/e  point  of  r^(x,0)/P^  and  the  y  direction  de- 
correlation  distance  is  given  by  the  1/e  point  of  T ^(0,y)/p These 
quantities  are 

Wlo  *  (3-41) 

and 

V'o  =  S'‘  Ls!GUxsi"2x*Gvxcos2xi'1/2  •  <3-42> 

Recall  that  under  the  frozen-in  approximation  the  drift  direc¬ 
tion  is  chosen  to  be  along  the  x  axis  and  the  decorrelation  time  of  the 
signal  is  related  to  the  decorrelation  distance  by  the  effective  velocity 
(Equation  2-87).  The  decorrelation  time  of  the  signal  out  of  the 
antenna  will  then  be  given  in  terms  of  the  decorrel ation  time  t0  of  the 
incident  signal  by 

ta/t q  *  *Ax/iQ  .  (Frozen-in  Approximation)  (3-43) 

The  filtering  of  the  decorrelation  time  is  a  reflection  of  the  fact  that 
signal  energy  arriving  at  large  angles  varies  more  rapidly  than  energy 
arriving  at  small  angles.  However,  under  the  turbulent  approximation,  the 
decorrelation  time  of  the  signal  is  independent  of  its  angle-of-arri val 
(see  Equation  2-100).  For  this  approximation,  there  is  no  filtering  of 
the  decorrelation  time  by  the  antenna  and 

t^/tq  =  1  .  (Turbulent  Approximation)  (3-44) 
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3.2.6 


Pouter  Impulse  Response  Function. 


The  power  impulse  response  function  gives  the  mean  power 
received  in  the  delay  interval  t  to  t  +  dx  from  a  transmitted  impulse.  At 
the  output  of  an  antenna,  this  function  is 

Ga(t)  =  (2*)"2  //  d2^  Gr(^)S(^,x)  .  (3-45) 

This  integral  cannot  be  done  in  closed  form  for  general  anisotropic 
scattering  and  antennas  and  for  arbitrary  values  of  the  parameter  a. 
However,  using  the  fact  that 

Limit  a  exp(-a2x2)  =  5p(x)  ,  (3-46) 

a  +  oe 

allows  an  analytic  result  to  be  obtained.  In  this  limit,  the  GPSD  of  the 
electric  field  incident  on  the  plane  of  the  receiver  becomes 

S(Kx,Ky,t)  =  (»//?)  exp[-(K*+a-2IC*)t*/4]  » 

sFf[“coh'-#(K(*Ky)eo/4'//?  1  '  (3-47) 

This  geometric  optics  limit  then  results  in  a  delta  function  relationship 
between  angle  and  delay  and  requires  that  x  be  greater  than  or  equal  to 
zero  in  order  for  the  GPSD  to  be  non-zero. 

The  power  impulse  response  function  for  general  anisotropic 
scattering  and  antennas  now  becomes 

(  ("coh/M)  exp(-g‘"cohT)  Io(92“coh')>  '  -° 
e  It\  =  '  (3-48) 

0  ,  T  <  0 
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where  I0  is  the  modified  Bessel  1  function  of  order  0  and  where 

91  =  C(62GUX+Gvy)c0s2x  +  (62GVX+Guy)sin2x]/(262A) 

92  *  [(6uy-626vx)2sin2x  ♦  (Gvy-52Gux)2cos2x]1/2/(262A)  .  (3-49) 

It  is  easy,  if  not  somewhat  tedious  to  show  that 
00 

/  GA(T)dT  =  PA  (3-50) 

which  is  equal  to  the  mean  received  power  at  the  output  of  the  antenna. 

In  the  limit  that  both  the  scattering  and  the  antenna  are  iso¬ 
tropic,  the  expression  for  the  power  impulse  response  function  reduces  to 

(  “coh  ex')(-G“coht)  ’  T  i  0 
Ga(t)  =  (3-51) 

(  0  ,  T  <  0 

where,  in  terms  of  the  standard  deviation  of  the  angle-of-arri val 
fluctuations  oq  and  the  3  dB  beamwidth  60, 

G  ■  1.0  +  (81n2)(cj0/0o)2  .  (3-52) 

Under  these  conditions,  the  quantity  G  is  equal  to  the  scattering  loss. 

The  cumulative  delay  distribution  C(t)  of  signal  energy  may  be 

defined  as 

C(t)  =  /  [GA(t)/PA]dt  .  (3-53) 

This  function  gives  the  fraction  of  received  signal  energy  arriving  with  a 
delay  less  than  or  equal  to  t.  For  isotropic  scattering  and  for  an  iso¬ 
tropic  antenna,  the  cumulative  delay  distribution  is 

C(x )  =  1.0  -  exp(-GurohT)  .  (3-54) 
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Without  an  antenna,  the  value  of  C(t)  is  equal  to  0.8  for  w  =  1.61. 
With  an  antenna,  as  the  value  of  6  or  the  scattering  loss  increases,  the 
80  percent  point  on  the  distribution  occurs  for  smaller  values  of  “C0^T* 

3.3  COMPARISON  OF  UNIFORMLY  WEIGHTEO  CIRCULAR  AND 
GAUSSIAN  ANTENNA  FILTERING  EFFECTS. 

The  filtering  equations  presented  in  the  previous  sections  are 
for  antennas  with  a  Gaussian  beam  profiles  which  do  not  have  sidelobes. 

It  is  of  interest  to  evaluate  the  filtering  equations  for  implementable 
antenna  beam  profiles  in  order  to  assess  the  effects  of  sidelobes  on  the 
statistics  of  the  signal  out  of  an  antenna  and  to  assess  the  accuracy  of 
the  Gaussian  beam  profile  filtering  equations. 

To  facilitate  this  calculation,  uniformly  weighted  circular 
antennas  and  isotrooic  scattering  will  be  assumed.  With  these  assump¬ 
tions,  the  moments  of  delay  defined  in  Equation  3-28,  the  equation  for  the 
spatial  coherence  function  (Equation  3-39)  and  the  cumulative  distribution 
of  the  received  signal  energy  (Equation  3-53)  can  all  be  written  in  terms 
of  a  set  of  single  integrals  which  can  be  evaluated  numerically.  From 
these  integrals,  the  scattering  loss  and  the  ratios  f^/fQ  and 
can  be  calculated.  These  quantities  are  a  function  only  of  the  ratio 
D/£0  where  D  is  the  diameter  of  the  circular  antenna. 

Under  isotropic  scattering  conditions  and  for  the  isotropic 
antennas  considered  here,  the  filtering  equations  for  Gaussian  antennas 
reduce  to 


where  the  function  G  is  defined  in  Equation  3-52  in  terms  of  the  standard 
deviation  of  the  angle-of-arri val  fluctuations  and  the  3  dB  antenna  beam- 
width.  For  uniformly  weighted  circular  or  square  apertures, 

1.0  +  0.265(D/£q)2  (Circular  Antenna) 

,  (3-56) 

1.0  +  0.358(0/ tQ)  .  (Square  Antenna) 

Figures  3-3,  3-4,  and  35  show  the  scattering  loss,  decorrelation 
distance,  and  frequency  selective  bandwidth  respectively  at  the  outputs  of 
Gaussian  (solid  lines)  and  uniformly  weighted  (dots)  circular  antennas. 
Figure  3-6  shows  the  value  of  w  for  which  the  cumulative  delay  distri¬ 
bution  is  equal  to  0.8. 

As  can  be  seen  from  the  figures,  there  is  good  agreement  in  the 
scattering  losses  and  the  decorrelation  distances  between  the  Gaussian 
approximation  and  the  uniformly  weighted  aperture  antennas  over  the  range 
of  D/£0  from  0.1  to  100  shown  in  the  figures.  The  scatterina  loss  is 
primarily  determined  by  the  fraction  cf  signal  energy  that  is  scattered 
out  of  the  main  beam  of  the  antenna  (Dana  1981)  and  is  insensitive  to  the 
sidelobe  structure  for  antennas  with  sidelobes  that  are  less  than  or  equal 
to  those  of  uniformly  weighted  apertures.  Under  conditions  where  D  »  lQ, 
an  antenna  must  only  be  displaced  approximately  a  distance  0  for  the 
signal  out  of  the  antenna  to  decorrelate.  Thus  under  these  conditions, 
the  decorrelation  distance  depends  primarily  on  the  aperture  size  and  is 
insensitive  to  the  sidelobe  structure. 

The  frequency  selective  bandwidth  of  the  signals  out  of  uniform¬ 
ly  weighted  circular  antennas  agree  with  the  Gaussian  approximation  values 
only  for  D/i0  less  than  3.  For  values  of  this  ratio  greater  than  3, 
there  is  less  filtering  of  the  frequency  selective  bandwidth  with  the 
uniformly  weighted  aperture  with  sidelobes  than  there  is  with  the  Gaussian 
antenna  without  sidelobes.  This  can  be  understood  by  considering  Equation 
3-23  for  the  ueiay  (numerics  of  the  received  siqnal.  The  GPSO,  Sn(K. ,,K t)  , 
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Figure  3-3.  Scattering  loss  for  Gaussian  and  uniformly  weighted 
circular  antennas. 
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Figure  3-4.  Decorrelation  distance  for  Gaussian  and  uniformly  weighted 
circular  antennas. 
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in  this  equation  couples  delay  and  angle  by  requiring  that  the  signal  < 

arriving  at  long  delays  also  arrive  at  large  angles.  The  moments  of  delay  i 

<tn>  then  become  more  sensitive  to  the  details  of  at  large  angles,  and  i 

therefore  to  the  sidelobe  structure  of  the  antenna,  as  n  increases.  As  | 

D/£0  increases,  more  of  the  signal  energy  is  scattered  out  of  the  main 
beam  and  Into  the  sidelobes  of  the  antenna.  Although  this  energy  contri¬ 
butes  little  to  the  total  received  power,  it  does  increase  the  values  of  ! 

the  delay  moments  of  the  output  signal  relative  to  those  of  the  signal  out  j 

of  an  antenna  without  sidelobes  and  thereby  decreases  the  frequency  selec-  1 

tive  bandwidth. 

i 

The  frequency  selective  bandwidth  is  intended  to  reflect  how  the  j 

bulk  of  the  signal  energy  arrives  in  delay.  To  see  how  the  uniformly 
weiohted  circular  antenna  affects  the  delay  distribution  of  the  received 
signal  energy,  the  80  percent  point  of  the  cumulative  delay  distribution 
(i.e.  the  value  of  U)C0^T  f°r  which  C(uco^t)  =  0.8)  is  plotted  in  Fiqure  3- 
6  for  both  Gaussian  and  uniformly  weighted  antennas.  Over  the  entire 
range  of  D/i0  values  from  0.1  to  100,  there  is  close  anreement  between 
the  Gaussian  antenna  values  and  the  uniformly  weighted  antenna  values 
indicating  that  at  least  80  percent  of  the  signal  energy  is  received  with 
a  similar  distribution  in  delay.  Even  for  the  largest  value  of  D/i0, 
more  than  80  percent  of  the  received  signal  energy  is  received  in  the  main 
lobe.  Thus  the  agreement  in  the  80  percent  values  is  a  reflection  of  the 
fact  that  a  Gaussian  function  provides  a  good  approximation  to  4 J i 2 (O/C 
within  its  main  lobe.  The  large  difference  in  the  filtered  values  of  the 
frequency  selective  bandwidth  ratio  seen  in  Figure  3-5  is  then  attributed 
to  large  differences  in  the  delay  distribution  of  no  more  than  20  percent 
of  the  received  signal  energy.  The  Gaussian  values  of  f^/f  will  there¬ 
fore  provide  an  accurate  description  of  the  delay  distribution  of  the  bulk 
of  the  received  signal  energy  and  should  be  used  for  a  1 1  antennas  provided 
a  reasonable  fit  can  be  made  to  the  main  lobe  usina  the  Gaussian  approxi¬ 
mation. 
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igure  3-6.  Cumulative  delay  distribution  80  percent  value  for  Gauss 
and  uniformly  weighted  circular  antennas. 


SECTION  4 

TRANSPONDER  COMMUNICATION  LINKS 


The  filtering  equations  for  a  transponder  communication  link 
with  four  antennas  and  two  uncorrelated  propagation  paths  will  be  pre¬ 
sented  in  this  section.  The  mean  power,  decorrelation  distances,  decor¬ 
relation  time,  and  frequency  selective  bandwidth  of  the  received  signal 
will  be  calculated.  In  order  to  simplify  the  calculations,  it  will  be 
assumed  that  all  antennas  are  isotropic  and  are  pointed  without  error, 
that  the  transponder  is  ideal  (i.e.  the  transmitted  signal  is  identical  to 
the  received  signal),  and  that  both  propagation  paths  suffer  strong  scat¬ 
tering  effects.  With  this  latter  assumption,  the  mutual  coherence  func¬ 
tion  derived  in  Section  2  applies  separately  to  each  path. 

The  bistatic  transponder  propagation  geometry  is  shown  in 
Figure  4-1.  Although  they  are  not  indicated  in  the  figure,  it  is  assumed 
that  both  the  uplink  from  the  sender  to  the  transponder  and  the  downlink 
from  the  transponder  to  the  receiver  have  transmitting  and  receiving 
antennas. 

The  mutual  coherence  function  of  the  signal  out  of  the  receivers 
antenna  wll1  be  used  to  calculate  the  filtering  equations.  With  the 
assumption  that  the  scattering  effects  on  the  two  propagation  paths  are 
uncorrelated,  the  mutual  coherence  function  can  be  written  down  directly 
as  the  product  of  the  coherence  functions  of  the  two  paths.  However,  some 
insight  can  be  gained  by  constructing  the  received  voltage  and  then  by 
calculating  its  mutual  coherence  function  so  this  procedure  will  be  used. 
The  frozen-in  approximation  will  be  used  to  model  the  temporal  fluctua¬ 
tions  and  to  calculate  the  decorrel ation  time  of  the  received  signal. 
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Figure  4-1.  Bistatic  transponder  link  geometry. 

4.1  MUTUAL  COHERENCE  FUNCTION. 

Consider  a  monochromatic  wave  Uu(C,w,t)  with  RF  angular  fre¬ 
quency  o)  that  has  propagated  along  the  uplink  path  from  the  sender  and 
that  Is  incident  on  the  plane  of  the  transponder  at  time  t  and  at  position 
t.  Tha  effects  of  the  sender's  transmitting  antenna  are  implicitly  con¬ 
tained  in  Uu.  The  voltage  at  the  output  of  the  transponder  receiving 
antenna,  without  pointing  errors,  can  be  written  as 

00 

eu(C,»,t)  -  //  Aru(i-!')Uu(!\<o,t)d2£’  (4-1) 

•  00 

where  Aru  is  the  aperture  weighting  function  of  the  receiving  antennna 
on  the  uplink.  If  a  monochromatic  wave  is  transmitted  by  the  transponder, 
then  the  incident  voltage  on  the  plane  of  the  receiver  at  position  c  after 
propagating  along  the  downlink  path  from  the  transponder  to  the  receiver 


69 


will  be  Ud(?>w,t).  Again,  the  effects  of  the  transponders  transmitting 
antenna  are  implicitly  included  in  Ud .  If,  however,  the  voltage  e(l,ui,t) 
is  transmitted,  then  the  incident  voltage  at  the  receiver  will  be 
Ud (c ,w,t)eu (l,w,t-tp)  where  t  is  the  nominal  propagation  time  from  the 
transponder  to  the  receiver.  The  voltage  at  the  output  of  the  receiving 
antenna  is 

ao 

ed(c,»,t)  =  H  Ard(c-c')Ud(c,,o.,t)eu(c,u,t-tp)  d2c'  (4-2) 

«*od 

where  Ard  is  the  aperture  weighting  function  of  the  receiving  antenna  on 
the  downlink. 

The  two-position,  two-frequency,  two-time  mutual  coherence  func¬ 
tion  of  the  received  signal  is  then  given  by  calculating  the  expectation 

rb(p»<j,d»tj)  =  <ed(ci  ,<D1,t1)ed*(C2»aj2tt2)>  (4-3) 

which,  for  statistically  stationary  processes,  must  be  a  function  only  of 
the  difference  d  s  q  .  C2»  “d  =  wl  "  "a*  and  td  =  ti  -  t2,  The  subscript 
b  denotes  the  coherence  function  of  the  bistatic  signal  that  has  propa¬ 
gated  from  the  sender  up  to  the  transponder  and  then  down  to  the  receiver. 

The  transponder  acts  as  a  point  source  for  the  receiver  just  as 
the  sender  is  a  point  source  of  the  signal  incident  at  the  transponder. 
Then  because  the  coherence  function  for  either  one-way  path  is  independent 
of  the  point  of  origin,  the  coherence  function  of  the  bistatic  siqnal 
depends  only  on  the  distances  z.  ( i -1 ,4)  and  on  p  and  is  independent  of 
the  positions  of  the  sender  or  the  transponder  as  long  as  the  z\  dis¬ 
tances  are  held  constant. 

At  this  point,  it  is  convenient  to  include  the  frozen-in  approx¬ 
imation  explicitly  in  the  formulation  by  writing  the  voltages  and  Ud  as 
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and 


uu(c.«,i)  =  uu(c-vut,^,o) 


(4-4a) 


ud(;,u.,t)  *  ud(c-vdt,w,0)  (4-4b) 

where  v  and  v  ,  are  che  effective  drift  velocities  of  the  random  diffrac- 

U  0 

tion  patterns  as  seen  respectively  at  the  transponder  on  the  uplink  and  at 
the  receiver  on  the  downlink. 


by  combining  Equotions  4-1  through  4-4,  the  coherence  function 
can  be  written  as 

OO  00 

rb(P.“d»td)  »  //d2^  A^U-^)  //  A * 

_  OD  _  OO 

<»u^i-vu(tl-tp),»1,o]uJ[C2---u(t2-tD),u.2,n]>  - 

oo  ao 

//  d2C|  ArdvCj-t;)  //  d2C'  Ard(C2-^)  * 


<ud(c;  vdt1,^1,0)Ud(q-vdt2,u2,0)>  .  (4-5) 

The  two  terms  in  the  brackets  <  >  are  recognized  as  the  £ync- 

tions  of  U  and  : '  . 

U  J 

Consider  for  the  moment  the  received  si  nal  wi‘  ecCional 

antennas  at  the  sender,  transponder,  and  r-ceiver.  The  received  signal 

will  then  be  just  the  product  U  Ud  and  bee  >e  the  propagation  effects  on 

the  two  paths  are  assumed  to  be  uncorre  1 -.ted,  =  r  Td#  As  was  done  in 

Section  2,  the  functions  r  and  r,  C3n  be  written  as  the  product  of  a  free 

u  o 

space  term  r,d  and  a  scintillation  term  Tj.  If  the  bis1".  c  ee  space 
term  is  set  equal  to  r  T  the  bi static  sc  in’  laci"'  .will  be 
T.  ,  -  r  ,1',.  In  what  follows,  only  the  contr.u,’'  r  .  .  tne  received 


voltage  due  to  scintillation  will  be  considered  and  an  expression  for 
wHl  be  derived  in  terms  of  the  uplink  and  downlink  GPSDs  and  the  four 
an  tennas . 

Equation  4-5  can  be  split  into  an  uplink  factor  times  a 

downlink  factor  r.  where  the  downlink  factor  is 
Ad 

00  CO 

rAd( P»“d»td)  =  JJ  d2^i  Ardf «!-«!)  //  Anj(vfy  * 

.00  .  00 

rdi[(ci-cp-:dtd.«d.o]  •  (*-6) 

The  subscripts  Ad  and  Au  are  used  to  denote  the  coherence  functions  at  the 
outputs  of  the  receiving  antennas  on  the  uplink  and  downlink.  It  is 
understood  that  rnj  and  r.  represent  the  coherence  functions  of  the  scin- 
ti llation  part  of  the  siqnal.  This  eauation  can  be  further  reduced  by 
writing  the  aperture  weiahtinq  function  in  terms  of  its  Fourier  transform 
(see  Equation  2-114).  After  changing  the  order  of  integration,  rfl^ 
becomes 

ao  co 

rAd(p»U)d»td )  =  (2")' 4  //  d2Kj  Arcj  ( K  i )  //  d2K2  A*j(K2)  x 

.00  -to 


a* 

exp[ iK i - ( ^ ! ) ]  fj  d2?^  exp[-iK2* (c2-C^) ] 

.  CO 


X 


rdn  4 i*C2'vdtd’U)d*u^ 

and  after  a  change  of  variables  from  c'j  tc  c  =  c1  -  r.'2  •  vgtg. 


/  a  ~r  \ 

\h-/  ; 


)Z 
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rAd(P."d.td)  =  (2*)'4  //  Ard(^)  exp[  ii<1*(c1-vdtd)] 

.00 


If  d2K2  A*d(K2)  exp  (-ii<2*52)  //  d2^  exp[ic2* (^-K^ ]  x 

_  OD  .00 


//  d2c  exp  ( - iK 1 • c)  rdi(c,ud,0) 

.00 


(4-8) 


The  last  integral,  which  will  be  denoted  (i<i ,0^ ) ,  is  the  Fourier  trans¬ 
form  of  the  angular-delay  part  of  the  GPSD  for  the  downlink  signal.  The 
subscript  A  is  attached  because  S ^  includes  the  effects  of  the  trans¬ 
ponder's  transmitting  antenna.  The  integral  over  results  in  the  delta 
function  ( ( k2 -K 1 ) .  With  these  substitutions,  the  coherence  function 
for  the  downlink  is 


=  {Zv)~2  U  exp  [ii<H*(p-vHtH)]  (4-9) 


d  rd  d  Ad  d  d ' 


d  d 


where  G  is  the  beam  profile  of  the  receiving  antenna  on  the  downlink, 
ro 


In  a  similar  fashion,  the  coherence  function  r^u  can  be 
derived.  However,  it  can  also  be  written  down  directly  from  Equation  4-9 
with  the  following  two  observat  ior- .  First,  the  bistatic  mutual  coherence 
function  is  the  expectation  of  received  voltage  for  two  positions  of  the 
receiver  but  for  only  one  position  of  the  transponder.  Thus  the  uplink 
contribution  to  contains  only  carrier  frequency  and  time  differences. 
Second,  the  expectation  of  U  U*  that  appears  in  Equation  4-5  depends  only 
on  the  argument  differences  and  is  therefore  independent  of  the  nominal 
propagation  time  t  .  The  uplink  coherence  function  is  then  given  by 


r  A  ( Co ,  %  1 1 
Au v  d  d 


)  »  (2ir) 


_2 


ff  d2£  6  (i<  )SA  (K  ,U.)  exp  (■ 

1 1  u  ru  u  Au  u  d 


iK  • v  t.)  (4-10) 

u  u  d' 


where  G  is  the  qain  function  of  the  receivinq  antenna  on  the  uplink  and 
,  ru 

SAu  is  the  Fourier  transform  of  the  angular-delay  GPSO  for  the  uplink  in¬ 
cluding  the  effects  of  the  sender's  transmitting  antenna. 


The  effect  of  a  transmitting  antenna  is  to  filter,  as  a  function 
of  angle,  the  angular-delay  GPSO  of  the  siqnal  at  the  receiver.  Thus  the 
angular-delay  GPSO  of  the  incident  signal  on  the  plane  of  the  receiver 
including  the  effects  of  the  transmitting  antenna  pointed  at  the  receiver 
is 


St(K,t)  =  Gt(<K)S(K,t)  (4-11) 

where  3(K,t)  is  the  GPSO  for  an  omnidirectional  transmi tter,  Gt  is  the 
beam  profile  of  the  transmitting  antenna  and  *  =  zr/zt  is  a  scale  factor 
that  transforms  angles  at  the  transmitter  to  angles  at  the  receiver  (see 
Figure  2-3  and  Equation  2-90). 


The  mutual  coherence  function  of  the  signal  out  of  the  receiver 
for  a  transponder  communications  link  is  now  given  by 

rbl(p,“d’V  =  rAu^d,t:d' rAd^p,u,d’td^ 


(2")'2  If  d2K  G  (k  K  )G  (K  )S  (K  ,id.)  exp  (-iK-v  t  )  x 
v  ij  tu  u  u  ru  u  u  u  d  "  "  ri 


u  u  d ' 


:2")'2  II  d2Kd  Ctd(«/d)Grd(Kd)Sd(Kd..d)  exp  ((iiV<Htd)] 


where  G^u  is  the  Deam  profile  of  the  transmitting  antenna  on  the  uplink 
and  Gtd  is  the  beam  profile  of  the  transmitting  antenna  on  the  downlink. 
The  <  scale  factors  are 


<U  =  22/Zl 


(4-13a) 


<d  *  24/z3  .  ( 4- 13b ) 

For  satellite  transponders,  the  distance  from  the  satellite  to  the  iono¬ 
sphere  will  usually  be  much  greater  than  the  distance  from  the  ionosphere 
to  the  ground.  Therefore,  >  1  and  <  1  under  most  circumstances. 

4.2  SCATTERING  LOSS. 

The  scattering  loss  suffered  on  the  transponder  communications 
link  depends  on  the  GPSDs  of  the  two  propagation  paths  and  on  the  beam 
profiles  of  the  four  antennas.  Under  the  assumption  that  the  two 
propagation  paths  are  uncorrelated,  the  total  scattering  loss  will  be  the 
product  of  an  uplink  contribution  and  a  downlink  contribution. 

The  mean  received  power  for  the  bistatic  path  is 

Pb  -  rbl(5,o.o)  *  PuPd 

00 

»  (2*)-2  II  d2Ki  6tu(Ku^l)Gru(itl)Su(it1,0)  x 

«oo 

ao 

( 2* )" 2  II  d2i<2  Gtd(<dK2)Grd(K2)Sd(K2,0)  (4-14) 
.00 

where  Py  is  equal  to  the  first  integral  and  where  P^  is  equal  to  the 
second  integral.  The  function  S^(K,0)  is 


Sc(K,0)  =  I  S?(K,T)dx 


■  <"VV  exp  l-(Kx'o?/4)  -  (KyV4S5>) 


(4-16) 
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where  S^(K,t)  is  the  angular-delay  part  of  the  GPSD  for  the  uplink  (?=u) 
or  the  downlink  (C=d).  In  the  developments  that  follow,  each  path  will 
have  its  own  decorrelation  distance  at  the  receiver  frequency 

selective  bandwidth  fQ?,  and  penetration  angle  and  scale  size  ratio  65. 
Also,  there  will  be  no  assumed  relationship  between  the  orientations  of 
the  x  axes  on  the  two  paths. 

The  antenna  beam  profiles  for  this  calculation  are  assumed  to  be 
isotropic  and  Gaussian  and  to  be  pointed  without  error.  The  combined 
transmitting  and  receiving  beam  profiles  may  then  be  written  as 

6tr(<i)G  _(i<)  =  exp  [-(«£< 2  +  «2  )K2]  exp  [-(«2  *2  +  aL)K2J  (4-16) 


JtO  C,vrV 


ri¬ 


te  C 


where  a  and  a  are  given  in  terms  of  the  transmitting  and  receiving 

Is  « 

antenna  beamwidths  respectively  by  Equation  3-3.  The  combined  beam 
profile  is  again  Gaussian  with  a  beamwidth  that  is  smaller  than  either  of 
the  two  antennas. 


The  scattering  loss  the  €  path  is 


LS?  ■  Pi‘  ‘  [Gx6Gy«l1/2 

(4-17) 

where 

r’x?  *  1  +  *  4ar5/4o6 

( 4- 13a ) 

and 

V  *  1  *  ■ 

(4- 18b) 

If  the  antennas  are  uniformly  weighted  circular  apertures,  then  the  G 
functions  become 

6xe  *  1  ♦  0-26S  («f<  VV'2  *  <V‘oC!21  (4-19a) 

and 

GyC  =  1  *  0.265  [<2(  VV2  *  ! W)2]  (',-19b) 
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where  and  Or^  are  the  transmitting  and  receiving  antenna  diameters 
respectively.  The  total  loss  of  the  bistatic  transponder  link  is 

LS  ■  LSuLSd  *  (GxuV\dGydT/2  ■  (4-20) 

The  decorrelation  distance  *0g  In  these  expressions  Is  that  of 
the  signal  incident  on  the  plane  of  the  receiving  antenna  when  an  omni¬ 
directional  signal  is  transmitted.  Recall  from  Section  2  that  the  decor¬ 
relation  distance  is  not  a  reciprocal  quantity.  If  an  omnidirectional 
signal  were  transmitted  at  the  receiver's  location,  the  decorrelation 
distance  of  the  signal  at  the  plane  of  the  transmitter  would  be  *0r/Kr* 

Cn  the  uplink  to  a  geosynchronous  satellite  transponder,  *  >  1  and  the 
uplink  scattering  loss  is  determined  primarily  by  the  loss  due  to  the 
transmitting  antenna.  However  on  the  downlink,  <d  «  1  and  the  downlink 
scattering  loss  is  determined  primarily  by  the  loss  due  to  the  receiving 
antenna.  In  either  case,  it  is  the  antenna  closest  to  the  ionosphere  that 
determines  the  scattering  loss  of  the  link. 

4.3  FREQUENCY  SELECTIVE  BANDWIDTH 

The  time  delay  jitter  of  the  bistatic  signal  at  the  output  of 
the  receiver  can  be  calculated  from  the  power  impulse  response  function 

co 

Pb(t)  =  (2*)-1  /  Tbl(5,u>d,0)  exp( ia>d-r )do>d  (4-21) 

a  00 

where  Ph(i)dT  is  the  fraction  of  signal  energy  arriving  in  the  delay 

D  f*  h 

interval  x  to  t  +  dx.  The  n  delay  moment  of  the  signal  energy  is 

00 

P  <Tn>  =  /  Tn  P  (x)dx  (4-22) 

where  P.  is  the  n  =  0  moment  and  is  equal  to  the  mean  received  power.  Th« 
o 

evaluation  of  the  delay  moments  is  simplified  by  considering  the  form  of 
P(x) .  The  power  impulse  response  function  can  be  written  as 


(4-23) 


Pfa(T)  =  (2*)"1  /  rA  (uid,0)r.d(5,u.,0)  exp(iajdT)d(Dd 

.40  ) 

This  integral  can  also  be  written  in  the  delay  domain  as  a  convolution: 

00 

Pb(T)  a  /  Pu(T‘)Pd(t-T')dT*  (4-24) 

.  OP 

where  P  (t)  and  Pd(T)  are  the  power  impulse  response  functions  of  the  up¬ 
link  and  downlink  paths  respectively.  These  functions  are 

<30 

Py(r)  *  ( 2 it ) ~ 1  /  rAu(^d,0)  exp( iaidT)d0Jd  (4-25) 

•OP 

and 

CD 

pd(*)  =  (2*)“l  /  rAd(^.ud.°)  exp(io>dT)dwd  .  (4-26) 

aQO 

Now  the  delay  moments  of  the  received  signal  can  be  evaluated  in 
terms  of  the  delay  moments  of  each  of  the  paths: 

00 

Pc<Tcn>  =  /  x n  Pc(x)dx  (4-27) 

«  00 

for  £  =  u  or  K  -  d.  The  mean  received  power  is 
00  00  00 

Pb  =  /  Pb(t)dT  =  /  dT  /  dT(P  (T')Pd(T-f)  (4-28) 

•  <3P  •  OO  .  00 

which  after  some  manipulations  reduces  to  the  previous  result,  given  in 
Section  4.2,  that  Pb  =  P^.  Using  this  formalism,  it  is  easy  to  show 
that  the  first  two  delay  moments  of  the  bistatic  signal  are 

<y  ■  <v  *  <v  (4-29) 


and 


<t2>  =  <t2>  +  2<t  XT.)  +  <t2> 


(4-30) 
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The  time  delay  jitter  of  the  bistatic  signal  is  then 


2  2  2 

a *  =  a  +  a  , 

Tb  tu  rd 


and  the  i  equency  selective  bandwidth  of  the  bistatic  signal  out  of  the 
receiver's  antenna  is 


fAb  ■  l/(2’°Tb>  • 


(4-32) 


The  problem  has  now  been  reduced  to  that  of  finding  the  fre¬ 
quency  selective  bandwidths  of  the  signals  at  the  outputs  of  the  receiving 
antennas  for  the  one-way  paths.  This  result  was  presented  in  Section 
3.2.3  for  one-way  paths  with  a  receiving  antenna  only.  However,  it  is  a 
simple  matter  to  generalize  those  results  to  include  the  effects  of  a 
Gaussian  transmitting  antenna.  The  frequency  selective  bandwidth  of  the 
signal  out  of  the  receiving  antenna  for  either  one-way  path  is 


fAC/fo« 


i <54G2  +  G2  l1/2 

1  Vxe  JySJ 


(4-33) 


where  5  *  u  or  5  «  d.  This  expression  is  valid  only  for  isotropic 
antennas.  The  quantity  f0^  is  the  frequency  selective  bandwidth  of  the 
signal  incident  on  the  plane  of  the  receiving  antenna  for  a  transmitted 
impulse  and  for  an  omnidirectional  transmitter  antenna.  The  time  delay 
jitter  for  either  path  is  given  by 


v  ■  I/(2'V 


(4-34) 


Finally,  the  frequency  selective  bandwidth  of  the  bistatic 
signal  at  the  output  of  the  receiving  antenna  is 


f  =  [f-2  +  f -  2  j - 1 / 2 
1  Ah  <■  '  All  1  AH  j 


(4-35) 


& 

if1 


ve  the  same  filtered  frequency 
1  to  fA^//2.  If,  however,  one  of 
selective  bandwidth  than  the  other 
of  fA?. 

AND  TIME. 

of  the  signal  at  the  receiver 
re  given  by  the  1/e  points  of 
ion  mutual  coherence  function  is 
iqual  to 

exp[-(«dy/*od)2/Gyd]  .  (4-36) 
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where 


t.  .  =  l  ,  /G  ./v  . 
Ad  od  xd  d 


(4-40) 


is  the  antenna  filtered  decorrelation  time  due  to  the  downlink  only  and 
where  v^  is  the  effective  downlink  drift  velocity  along  the  downlink  x 
axis.  A  similar  set  of  expressions  hold  for  the  uplink  two-time  coherence 
function  and  decorrelation  time: 


where 


rAu(5*°’td)  =  Pu  exP[-<vutd)2/(6xut5u)J 
=  Pu  exPt’^d/rAu)2^ 


(4-41) 


t  -l  /G  /v 
Au  ou  xu  u 


(4-42) 


is  the  antenna  filtered  decorrelation  time  due  to  the  uplink  only  and 
where  vu  is  the  effective  uplink  drift  velocity  along  the  uplink  x 
axis.  The  bistatic  path  has  a  two-time  mutual  coherence  function  given  by 


rAb(^,0*td>  =  rAu^*°*td^rAd^*0,td^ 

*  Vd  ex»i-<Ts* +  *Ad)‘3i  •  <4-43> 


The  bistatic  path  decorrelation  time  under  the  frozen-in 
approximation  is  then  given  by 


T 


Ab 


(Frozen-in  Approximation) 


(4-44) 


The  decorrelation  time  for  the  turbulent  approximation  can  be  obtained 
from  v.his  expression  by  noting  that  under  this  approximation,  the 
decorrelat;on  time  is  not  filtered  by  the  antennas.  Therefore,  the 
bistatic  decorrelation  time  becomes 


Tob  =  ^Tou  +  Tod^"l//2  ’  (Turbulent  Approximation)  (4-45) 
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where  t  and  t  ,  are  the  uplink  and  downlink  decorrelation  times  without 
ou  od 

antenna  effects.  This  expression  is  valid  only  if  the  Doppler  frequency 
spectrum  is  Gaussian.  If  this  is  not  the  case,  then  a  slightly  different 
relationship  holds  between  the  bistatic  path  decorrelation  time  and  the 
decorrelation  times  of  the  two  Individual  paths  (Dana  1982). 

If  both  the  uplink  and  the  downlink  have  the  same  decorrelation 
time  at  the  outputs  of  the  receiving  antennas,  then  under  the  frozen-in 
(turbulent)  approximation,  is  equal  To  f^A'2  (T0£/'/2)-  ^ 

however,  one  of  the  links  has  a  much  smaller  decorrelation  time  than  the 
other,  then  t  ( rQb)  is  equal  to  the  smaller  value  of  TA?(’0?)- 
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SECTION  5 

CHANNEL  SIMULATION 


The  purpose  of  this  section  is  to  describe  a  statistical  channel 
simulation  technique  that  allows  realizations  of  the  impulse  response 
functions  to  be  generated  at  the  outputs  of  multiple  antennas  with  spatial 
and  temporal  correlation  properties  given  by  the  6PS0  derived  in  Section  2 
and  with  Rayleigh  amplitude  statistics.  The  realizations  generated  with 
this  technique  represent  only  the  diffractive  part  of  the  received  voltage 
and  they  are  valid  only  under  strong  scattering  conditions  where  the  GPSD 
is  valid  and  where  Rayleigh  statistics  apply.  Under  these  conditions 
however,  they  represent  a  solution  of  Maxwell's  equations  for  propagation 
of  RF  waves  through  randomly  structured  ionization. 

The  basic  formalism  to  generate  statistical  realizations  of  the 
channel  impulse  response  function  without  antenna  effects  explicitly  in¬ 
cluded  was  developed  by  Wittwer  (1980)  for  isotropic  irregularities  and 
was  extended  by  Knepp  (1982)  to  the  case  of  elongated  irregul arities . 

(The  elongated  case  corresponds  to  a  90°  penetration  angle  and  to  an 
Infinite  axial  ratio.)  The  channel  simulation  technique  has  been  general¬ 
ized  in  this  report  to  include  the  effects  of  anisotropic  scattering  and 
the  effects  of  multiple  receiving  antennas  and  has  resulted  in  a  new  ONA 
software  channel  simulator  for  antenna  applications  (Dana  and  Wittwer, 
1985). 
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It  will  be  assumed  in  this  section  for  the  sake  of  limiting  the 
discussion  that  the  multiple  antennas  are  Identical  In  beam  profile,  ori¬ 
entation,  and  pointing  angle.  It  is  also  assumed  that  the  antenna  centers 
are  col  inear.  These  restrictions  may  be  easily  lifted  by  using  the  more 
general  formalism  reported  by  Wittwer  (1986). 

The  frozen-in  approximation  will  be  used  in  most  of  the  develop¬ 
ments  in  this  section  to  describe  the  temporal  fluctuations  of  the  inci¬ 
dent  signal.  The  generalization  of  the  channel  simulation  technique  to 
the  turbulent  approximation  is  discussed  in  Section  5.2  along  with  a 
second  technique  proposed  by  Wittwer  (1986)  that  is  valid  in  the  a  equal 
infinity  limit.  A  comparison  of  realizations  generated  under  the  frozen- 
in  and  turbulent  approximations  is  made  in  Section  5.4. 

The  scattering  of  the  signal  is  described  in  the  x-y  coordinate 
system  shown  in  Figure  5-1,  which  is  normal  to  the  1 ine-of-sight.  The  x 
axis  corresponds  to  the  direction  with  the  minimum  decorrelation  distance 


Antenna 

Center 


Figure  5-1.  Scattering  and  antenna  coordinate  systems. 
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which,  to  be  conservative,  is  chosen  to  be  the  direction  of  the  effective 
velocity  between  the  antennas  and  the  random  diffraction  pattern.  The 
antennas  are  described  in  the  u-v  coordinate  system,  which  lies  in  the 
x-y  plane,  where  the  antenna  centers  lie  along  the  u  axis.  The  angle 
x  (0  £  x  <  tt/2 )  is  the  rotation  angle  between  the  scattering  coordinate 
system  and  the  antenna  coordinate  system.  In  the  figure,  one  of  the 
antennas  is  depicted  at  the  origin  of  the  x-y  coordinate  system.  The 
antenna  is  specified  by  Pm  where  Pm  is  a  position  vector  in  the  x-y 
plane  that  points  to  the  mth  antenna  center. 

5.1  GENERATION  OF  REALIZATIONS  (FROZEN- IN  APPROXIMATION). 


The  technique  for  generating  the  realizations  of  the  impulse 
response  function  at  the  outputs  of  multiple  antennas  is  outlined  below. 
For  a  given  signal  delay,  random  samples  of  the  impulse  response  function 
at  the  outputs  of  the  antennas  are  generated  on  the  two  dimensional  x-y 
grid  shown  schematically  in  Figure  5-2.  The  y  coordinates  of  the  grid 
correspond  to  the  y  coordinates  of  the  multiple  antenna  centers  ym.  The 
coordinates  of  the  grid  correspond  to  the  x  positions  of  the  antenna 
centers  x„  as  a  function  of  time  due  to  the  effective  velocity  of  the 
antenna  relative  to  the  random  diffraction  pattern. 


5.1.1  Discrete  Evaluation  of  the  GPSD. 

The  first  step  in  generating  the  impulse  response  function  at 

the  output  of  an  antenna  is  the  evaluation  of  the  GPSD  on  a  discrete 

K  -  K  -  t  grid.  The  angular-delay  grid  centers  are  defined  by 
x  y 

Kx  ■  k4Kx 

Ky  = 


Antenna  Center 


where  k,  1,  and  j  are  integers.  The  starting  delay  ts  and  the  angular 
and  delay  grid  sizes  will  be  defined  in  Section  5.1.3.  The  mean  signal 
energy  at  each  grid  point  is 

(k+l/2)AKx  { 1+1/2) AKy  Tj+At/2 

Ek,U  *  {U)~Z  i  dKx  /  dKy  /  dT  SA(Kx*Ky’T)  (5"2) 

( k-1/2) AKX  (1-1/2)AK  t^-At/2 

where  S^(Kx,Kw,t)  is  given  by  Equation  3-22  in  the  x-y  coordinate  system 
The  mean  signal  power  in  the  jtf1  delay  bin  is 


and  the  total  mean  power  in  the  grid  is 

V^j  •  ‘5-4> 

j 

These  two  quantities  are  useful  in  verifying  the  statistics  of  random 
realizations  as  discussed  in  Section  5.3. 


oc 


L»l»] 


The  integral  over  delay  in  Equation  5-2  can  be  done  in  closed 
form  with  the  result  that  the  mean  signal  energy  in  the  k-l-j  grid  cell  is 

(k+l/2)AKx 

Ek,l,j  =  /  dKx  exp  (-K*£2/4)  x 

(k-l/2)AKx 

(l+l/2)AKy 

!  dKy  exp  (-kV/462)  G(Kx,Ky)  Ej(Kx,Ky)  (5-5) 
(1-1/2)ak 


where 


:j(Kx,Ky)  *  erf  {o[“>c0h(y^/2)  -  A(K2+Kj)  t2/4] HZ\  - 


erf  lCtua*coh(Tj -At/2)  -  A(K2+Kj)i*/4]//2} 


(5-6) 


and  where  erf  is  the  error  function.  Because  the  antenna  beam  profile 
appears  explicitly  in  the  integrand  of  Equation  5-5,  £.  ,  .  is  valid  only 

K  f  I  ,  j 

for  an  antenna  with  beam  profile  G(K  ,K  ).  However,  at  this  point  the 

x  y 

location  of  the  antenna  is  arbitrary  as  long  as  the  beam  profile  remains 

fixed.  If  the  different  antennas  have  different  beam  profiles  (or  the 

same  profile  with  different  pointing  angles),  then  a  separate  E.  ,  .  must 

be  calculated  for  each  different  G(K  ,K  ). 

x  y 


The  value  of  a  In  Equation  5-6  Is  chosen  to  be  large  enough  so 
that  under  the  most  disturbed  propagation  conditions  the  antenna  filtering 
of  the  frequency  selective  bandwidth  is  not  limited  by  the  value  of  a. 

The  procedure  for  choosing  a  is  as  follows.  The  ratio  f^/f  under  con¬ 
ditions  where  the  decorrelation  distance  approaches  zero  is  obtained  by 
combining  Equations  3-35  and  3-36  with  the  result 


l imit  f ^/f  -  a2  +  l 


(5-7) 


♦ 


Using  the  minimum  decorrel ation  distance,  the  maximum  value  of  f./f  is 

A  0 

obtained  from  Equation  3-38  for  a  set  equal  to  infinity.  The  value  of  a 
is  then  chosen  so  that  the  difference  between  a2  +  1  and  the  maximum  value 
of  f^/f  is  small.  As  discussed  by  Knepp  (1982),  it  is  possible  to  set  a 
to  infinity.  However,  the  resulting  discontinuities  in  the  6PS0  present 
numerical  difficulties  in  the  signal  generation  process  that  are  avoided 
with  a  finite  value  of  a.  An  alternative  channel  simulation  technique 
which  depends  on  a  being  infinite  is  discussed  in  Section  5.2.2.  This 
latter  technique  was  developed  to  generate  signal  realizations  under  the 
turbulent  approximation  but  it  could  be  applied  to  the  generation  of 
frozen-in  approximation  realizations  as  well. 

5.1.2  Random  Realizations. 

The  impulse  response  functions  will  be  generated  one  delay  at  a 
time  starting  with  the  smallest  delay  (j=0)  and  working  to  the  largest 
value  of  delay  ( j  = j  ) .  For  each  delay,  the  steps  that  are  performed 

iMaX 

are  outlined  in  this  subsection. 


The  Fourier  transform  h^  of  the  impulse  response  function  at 
the  output  of  an  antenna  is 

1  (4’2/4K/Vt)5k,l./fMTj  •  (5-8) 

The  normalization  factor  ( 4-rr2/ AK  AK  At)  has  been  chosen  so  that  after  h. 

A  Jr  H 

is  Fourier  transformed  from  the  angular  K  -  K  domain  to  the  spatial 

x  y 

domain,  h^AT  will  represent  the  voltage  received  during  the  delay  interval 
Tj  -  Ax/2  to  Tj  +  At/2.  The  quantity  ^  .  is  a  complex,  zero  mean, 
Gaussian  random  variable  with  the  properties 

<5.  ,  .£*  .  >  =  6,  6.  „6  . 

k,  1  ,J  a,8,Y  k,a  1,6  J ,Y 


"k.u"  ”° 

<5k.l,A,ey  2  0 


(5-9) 


OO 


where  6  is  the  Kronecker  delta  symbol.  The  random  numbers  5  may  be 
qenerated  using  the  equation 

C  =  «Mn  uj  exp(2^iu2)  (5-10) 

where  ui  and  U2  are  independent  random  numbers  uniformly  distributed  on 
the  interval  [0,1). 

There  is  no  explicit  y  coordinate  since  realizations  will  be 
generated  only  at  y  coordinates  corresponding  to  antenna  centers.  This  is 
accomplished  by  first  generating  N  samples  of  h.  in  the  K  domain  and 

-  y  ^  y 

then  by  Fourier  transforming  h^  to  the  y  domain  using  discrete  Fourier 
transforms  (DFTs).  One  OFT  's  performed  at  each  coordinate  grid  point 
for  each  unique  antenna  center  y  coordinate.  In  continuous  notation,  this 
transform  is 


Vm  Kx-T>  s  <1/2*>  i  h’A<Kx'VT)  eXP  (iV»)4K> 


(5-11) 


where  ym  is  the  y  coordinate  of  the  mt!l  antenna  center.  The  OFT  equiva¬ 


lent  of  this  equation  is 


Ny/2-1 


hA,m(^Kx,Tj)  =  (&Ky/2ir)  l  exp  (iUKyym)  hA(kAKx,UKy,T ,•)  .  (5-12) 

1— N  /2 


Mow  the  Kx  to  x  Fourier  transform  is  performed  for  each 
antenna.  In  this  case,  a  fast  Fourier  transform  (FFT)  may  be  used.  In 
continuous  notation 


h n  _(x,T)  =  (1/2*)  /  exp  (i«„x)  ha  m(Kv,x)dK 


A  m'S' 


(5-13) 
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The  discrete  equivalent  of  this  equation  is 


V2'1 

ha  fJkAx’Ti)  =  l  exp  ( ik*  AK  kAx)  h.  (k'AK  ,t.)  (5-14) 

A.m  j  x  k,=.N)_/2  *  A.m  x  j 

where  Nx  is  the  number  of  x  coordinate  grid  points. 

The  quantity  h^  m(kAx,T^)  is  the  impulse  response  function  at 
the  output  of  the  mth  antenna  when  that  antenna  is  located  at  x  position 
kAx.  Under  the  frozen-in  approximation,  a  time  step  At  may  be  associated 
with  Ax  using  the  effective  velocity: 


At  =  Ax(tq/*o) 


(5-15) 


where  t  is  the  decorrelation  time  of  the  incident  signal.  The  quantity 
o 

h^  m(kAx,tj)  then  also  represents  the  impulse  response  function  at  the 
output  of  the  mth  antenna  at  time  kAt. 


5.1.3  Grid  Sizes. 


The  grid  sizes  of  the  channel  simulator  are  chosen  on  the  basis 
of  the  statistics  of  the  signal  out  of  the  antennas  which  are  given  by  the 
filtering  equations  in  '"action  3.2.  The  angular  grid  sizes  are  chosen  so 
the  antenna  filtered  GPSD  will  be  small  or  zero  at  the  edges  of  the  angu¬ 
lar  grid  and  so  there  will  be  a  sufficient  number  of  samples  per  decor¬ 
relation  distance  to  resolve  the  fades  of  the  Rayleigh  amplitude  distribu¬ 
tion.  The  delay  grid  size  at  is  chosen  on  the  basis  of  the  sampling 
period  of  the  receiver  and  the  number  of  delay  samples  of  the 
realizations  is  chosen  so  the  total  delay  spread  of  the  realization  NpAt 
will  encompass  at  least  95  percent  of  the  delayed  signal  energy. 


(5-16) 


5. 1.3.1  The  Ax  -AK^  Grid.  The  x  grid  spacing  is 
AX  =  AAx/N0 

where  i ^  is  the  filtered  x  direction  decorrelation  distance  (Equation  3- 
41).  The  quantity  N0  is  the  number  of  samples  per  decorrelation  dis¬ 
tance  or  time  where  N0  should  have  a  value  of  10  or  larger  in  order  for 
the  realization  to  accurately  represent  the  duration  and  separation  of  the 
fades  of  the  Rayleigh  amplitude  distribution  (Dana  1982). 

The  number  of  x  grid  points  Nx  is  subject  to  two  constraints: 
First,  Nx  should  be  a  power  of  2  so  a  fast  Fourier  transform  can  be  used 
to  transform  from  the  Kx  domain  to  the  x  domain.  Second,  Nx  should  be 
greater  than  or  equal  to  100  N0.  Equivalently,  there  should  be  at  least 
100  x  direction  decorrelation  distances  in  the  realization  in  orde<*  for 
the  realization  to  represent  a  reasonable  sample  of  the  random  process. 
Under  the  frozen-in  approximation,  100  decorrelation  distances  corresponds 
to  an  elapsed  time  of  100  decorrelation  times. 

Here  reasonable  is  defined  in  terms  of  the  application.  Typi¬ 
cally,  the  realizations  are  used  to  exercise  simulations  of  transiono- 
spheric  communications  links  and  each  realization  with  N  /N  -  100  can  be 

X  0 

used  to  calculate  receiver  performance  averaged  over  100  decorrelation 
times.  The  resulting  receiver  performance  measures  will  have  some  statis¬ 
tical  variation  due  to  the  finite  time  duration  of  the  simulation,  This 
variation  can  be  reduced  by  either  using  1 o  er  realizations  or  by  using 
several  realizations  with  the  same  signal  ,  .meters  but  with  different 
random  numbers  and  by  averaging  the  results  from  the  multiple  realiza¬ 
tions.  This  latter  approach  allows  an  estimate  of  the  statistical  varia¬ 
tion  in  the  receiver  performance  to  be  made. 
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Returning  to  the  grid  spacing,  the  x  direction  distance  spanned 
by  the  realization  is  lx  =  NyAx  and  A«x  is 


AKx  =  2it/l_x  =  2it/NxAx 


(5-17) 


The  functional  dependence  of  the  GPSD  on  Kx  is  given  auroximately  by 

c  -  expHK^/2)2]  .  (5-18) 

In  discrete  notation,  Kx  =  kAKx  where  |k|  has  a  maximum  value  of  Nx/2 . 
Using  the  maximum  value  for  Kx  of  NxAKy/2,  the  value  of  t  at  the  edge 
of  the  Kx  grid  is  exp[-(nNQ/2) 2]  which  is  essentially  zero  for  Nq  equal 
to  10. 


5. 1.3. 2  The  AKy  Grid.  The  number  of  Ky  grid  points  is  determined  by 
requirinq  that  the  GPSD  is  small  at  both  ends  of  the  K  grid.  The  func- 

y 

tional  dependence  of  the  GPSO  on  Ky  is  given  approximately  by 


e  «  exp[-(K  A.  /2) 21 


(5-19) 


A  convenient  minimum  value  of  e  is  exp(-n2)  =  5.2  *  10‘°.  Using  this 
value  results  in  the  following  condition  on  Ny: 


Ny  >  W(AKy*Ay) 


(5-20) 


It  is  desirable  to  minimize  the  number  of  Ky  samples  generated  to  mini¬ 
mize  computer  memory  and  execution  time  requirements  for  signal  genera¬ 
tion.  A  reasonable  minimum  value  for  Ny  is  32. 


The  AKy  grid  spacing  is  then 


AK  =  2u/l 

y  y 


(5-21) 


where  the  length  of  the  realization  in  the  y  direction  is  chosen  as 


L  =  maximum  [16*  ;  4  maximum  (|y  i\] 


(5-22) 


The  first  condition,  L  =  16 jL  ,  results  from  letting  N  =  32  in 

*  ^  .  y  . 

Equation  5-20.  The  second  condition,  =  4  maximum  ( j ym | ) ,  is  chosen 
to  minimize  aliasing  of  the  impulse  response  functions  of  the  antenna  with 

the  largest  value  of  y  and  the  antenna  with  the  smallest  value  of  y  .  If 

rn  m 

the  value  of  L  is  determined  by  this  latter  condition,  then  the  value  of 
Ny  will  need  to  be  larger  than  32  in  order  to  satisfy  the  condition  on 
given  by  Equation  5-20. 

5. 1.3. 3  The  Delay  Grid  At.  The  delay  grid  size  is  usually  chosen  so 
there  are  an  integer  number  of  delay  samples  in  a  sample  period  of  the 
communications  link  receiver.  There  should  be  at  least  two  delay  samples 
per  symbol  in  order  to  represent  the  frequency  spectrum  of  the  symbols. 

The  number  of  delay  samples  is  chosen  so  the  delay  spread  of  the  realiza¬ 
tion  NqAt  will  encompass  at  least  95  percent  of  the  delayed  signal 
energy.  However,  the  delay  spread  of  the  signal  depends  on  the  anisotropy 
of  the  signal,  on  the  value  of  a,  and  on  the  antenna  filterinq.  For  a 
given  frequency  selective  bandwidth  of  the  incident  signal,  the  delay 
spread  of  the  signal  out  of  an  antenna  increases  as  the  antenna  filtering 
is  reduced  (or  as  Z0  is  increased),  as  the  value  of  a  is  decreased,  and 
as  the  penetration  angle  is  reduced  (as  the  incident  signal  becomes  more 
isotropic).  Without  antenna  filtering,  95  percent  of  the  signal  energy 
under  isotropic  scattering  conditions  arrives  with  delays  in  the  range 
-0.25  <  w  <  3.45  when  a  is  set  at  its  minimum  value  of  4  (Wittwer 
1980).  To  be  conservative,  this  same  criterion  is  applied  under  all 
conditions  The  minimum  delay  spread  of  the  impulse  response  function 
realizations  must  therefore  satisfy  the  condition 

-0.25  <  2nfftT  <  3.45  .  (5-23) 

The  minimum  number  of  delay  samples  required  is  then 


Nd  >  1  +  3.7/(2nfAVc) 


(5-24) 


where  the  first  delay  bin  (j=0)  has  a  center  delay  of 

t$  =  -0.25/(2*fA)  .  (5-25) 

5.2  GENERATION  OF  REALIZATIONS  (TURBULENT  APPROXIMATION). 

The  GPSD  developed  in  Section  2  has  the  form  given  by  Equation 
2-102  which  is  reproduced  here: 

S(Kj_,T,a,0)  =  SD(«0)S(Kpx)  .  (2-102) 

Under  the  turbulent  approximation,  the  Doppler  frequency  spectrum  S^w^) 
is  independent  of  the  angular-delay  part  of  the  GPSD  S(K^,t).  The  Doppler 
frequency  spectrum  and  therefore  the  decorrelation  time  of  the  signal  is 
then  the  same  for  any  angle  and  any  delay. 

Two  techniques  for  generating  realizations  of  the  impulse 
response  function  under  the  turbulent  approximation  will  be  discussed  in 
this  section.  The  first  technique  is  a  generalization  of  the  technique 
used  for  frozen-in  approximation  realizations  and  is  valid  for  appropri¬ 
ately  chosen  large  but  finite  values  of  the  parameter  a.  The  second  tech¬ 
nique,  proposed  by  Wittwer  (1985),  is  valid  only  for  infinite  values  of 
a.  However,  for  reasons  described  in  Section  2.8.1,  this  limitation  does 
not  reduce  the  usefulness  of  Wittwer's  technique.  Indeed,  it  has  the 
advantage  of  requiring  less  computer  storage  and  cpu  time  to  generate 
realizations  than  are  required  with  the  first  technique.  Both  techniques 
are  somewhat  different  than  that  described  in  Section  5.1. 

5.2.1  Finite  <*  Technique. 

The  starting  point  of  this  generation  technique  is  Equation  5-3 
which  is  rewritten  as 
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h.  (kAK  ,1AK  ,T,,t  )  =  (4n2/AK  AK  At)  5.  ,  .(t  )/E.  ,  . 
A  x  y  J  n'  x  y  k ,  1 ,  j v  n'  k ,  1 ,  j 


(5-26) 


where  tn  are  discrete  times  for  which  the  impulse  response  functions 
will  be  calculated.  The  temporal  variation  of  the  impulse  response 
function  is  given  by  the  temporal  variation  of  the  random  numbers 

l  ,-(0*  These  random  numbers  must  still  satisfy  the  conditions  given 
by  Equation  5-9.  The  grid  sizes  AK^  and  At  are  calculated  using  the  cri¬ 
teria  given  in  Sections  5. 1.3. 2  and  5. 1.3. 3  respectively.  The  grid  size 
AK^  for  this  model  is  chosen  using  the  same  arguments  that  were  used  to 
select  AKy.  That  is,  the  number  of  «x  grid  points  is  first  chosen  to 
be  32  and  then  Equations  5-21  and  5-22  (with  y  replaced  by  x)  are  used  to 
calculate  ak  . 

A 

The  impulse  resoonse  functions  at  the  outputs  of  multiple 
antennas  are  calculated  in  a  manner  similar  to  that  outlined  in  Section 
5.1.2.  At  a  given  time,  the  impulse  response  function  for  the  mth 
antenna,  with  center  coordinates  x^  and  y  ,  is  calculate^  from  hA  using  a 
two  dimensional  Fourier  transform.  In  continuous  notation, 


-  <2’>-2  I  dKx  >  dKyexP  (tKxxm*iV«i)  V^X-V’  •  (5'27) 


.oo  .00 


The  discrete  equivalent  of  this  equation  is 

h.  (t  .,t  )  =  (AK  AK  /4  tt  2 )  X 
A, m'  j ’  n 7  v  x  y  ’ 


Nx/2-l  N  /2-1 

I  I  exp  <1k«xxm*11«  yj  hA(k«  ,1«  ,T j,tn) 

k=-N  /2  l«-N/2 

*  J 

This  operation  is  repeated  for  each  set  of  unique  antenna  center 
coordinates  and  then  for  each  delay. 


(5-23) 
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At  the  next  time  step  t  =  t  +  At  (where  At- t  /N  ),  the 

n+1  n  oo 

random  numbers  r.  are  recalculated  according  to  the  Doppler  frequency 
K  »  *  y  J 

spectrum.  A  convenient  form  for  the  Doppler  spectrum  is  an  f-4  Dower 

spectrum  because  it  can  be  readily  synthesized  using  two-pole  low-pass 

filters.  The  noise  sample  ,  -(t^.,)  can  then  be  obtained  directly 

x » I , j  n+i 

from  5.  ,  .(t)  and  a  complex  white  Gaussian  noise  sample.  Once  the  noise 
k  » i ,  j  n 

samples  r,  have  been  updated,  the  impulse  response  functions  at  time  tn+^ 
are  calculated  using  the  above  equations. 


5.2.2 


Infinite  a  Technique. 


This  technique  will  utilize  the  delta  function  relationship 
between  angle  and  delay  given  by  Equation  3-47  for  the  infinite  ot  limit. 
For  the  jth  delay  bin,  the  GPSD  will  be  non-zero  only  for  the  angles 
that  satisfy  the  condition 


W’j"*'2’  <  •'<Kx*l<y),o/4  <  “coh<Vit/2> 


(5-29) 


In  this  limit,  o)CGh  is  related  to  the  frequency  selective  bandwidth  as 


w  h  -  2nf 
coh  0 


(5-30) 


The  impulse  response  function  may  now  be  generated  in  the  angular  domain 
without  regard  to  delay.  Once  this  has  been  done.  Equation  5-29  may  be 
used  to  assign  an  annulus  in  the  K  -  K  plane  to  each  delay  bin. 

X  V 

The  starting  point  for  this  technique  is  then  the  energy  in  a 

K  -  K  grid  cel  1 : 
x  y  3 


;k+l/2)AK  ( 1+1A2)  AK 

*  y 


Ek,l  ■  ^  >  *,  -  y 

(k-l/zrac,  ( 1-1/2)  :i< 

*  y 


;  dK  G(K  ,K  )S(K  K  ) 
x  y  x  y 


(5-31) 
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where  the  angular  power  spectral  density  is 


S(Kx,Ky)  »  /  S(Kx,Ky,t)dT 

-.00 

*  (TtJ/6)  exp  [-(K2t2/4)  -  (Kjt20/462)]  (5-32) 

The  random  angular  spectrum  of  the  signal  at  discrete  time  tn  is 
rewritten  from  Equation  5-26  as 

hA(k&KxlAKy,tn)  =  (4n2/AKxAKyKkj(tn)/Ekj  .  (5-33) 

For  the  jth  delay  bin,  the  Fourier  transform  of  the  impulse  response 
function  a':  time  tn  is 

!h^(kAKx,UKy,tn)/&i  ,  if  Equation  5-29 

(5-34) 

0  ,  otherwise  . 

The  impulse  response  functions  at  the  outputs  of  multiple  antennas  are 

then  generated  from  h  (kA«  ,1AK ,T.,t  )  using  Equation  5-28. 

m  x  y  j  n 

The  savings  in  computer  resources  between  this  technique  and 

that  described  in  Section  5.2.1  is  the  result  of  the  fact  that  only  the 

random  angular  spectrum  instead  of  the  random  angular-delay  spectrum  is 

required  at  each  discrete  time  for  the  latter  technique.  Thus  E^  ^  is 

stored  rather  than  F  .  and  at  each  discrete  time  the  random  numbers 

K  *  '  *  J 

5.  ,(t  )  are  updated  and  stored  rather  than  S.  ,  .(t  ). 

For  this  technique,  the  mean  signal  power  in  the  grid  is 

PG  “  -  ^  Ek  ,  (5-35) 

k  1 
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and  the  mean  power  in  the  j th 


delay  bin  is 


T  ,+  At/2 

J 

/ 

Tj- At/2 


6A(x')dV 


(5-36) 


where  Ga(*)  Is  the  power  Impulse  response  function  given  in  Equation 
3-48.  For  general  anisotropic  scattering  and  antennas,  this  integral  must 
be  performed  numerically. 


5.3  REALIZATION  PARAMETERS. 


A  check  can  be  made  of  the  generated  realizations  of  the  Impulse 
response  function  by  computing  the  mean  signal  power  in  the  grid,  scatter¬ 
ing  loss,  frequency  selective  bandwidth,  x  direction  decorrelation  dis¬ 
tance,  and  the  cross  correlation  between  antennas  and  then  by  comparinq 
these  quantities  with  their  ensemble  values  wnicn  are  qiven  In  Section 
3.2.  The  realization  signal  power  in  each  delay  bln  may  also  be  compared 
with  the  ensemble  values  to  ensure  that  the  realization  has  the  proper 
delay  distribution  of  signal  energy.  For  reasons  discussed  below,  this  is 
often  a  better  check  than  comparing  the  realization  and  ensemble  frequency 
selective  bandwidths. 


Because  each  realization  contains  a  finite  number  of  decorrela¬ 
tion  distances  or  times,  the  agreement  between  the  realization  parameters 
and  the  ensemble  values  will  not  be  exact.  In  addition,  the  realization 
parameters  will  also  vary  from  antenna  to  antenna.  However,  as  the  reali¬ 
zations  are  made  longer  in  distance  (i.e.  Nx  is  increased  while  NQ  is  kept 
fixed)  the  agreement  will  improve. 

An  exception  to  this  is  the  frequency  selective  bandwidth  which 
requires  both  a  long  realization  in  distance  or  time  and  a  fine  gridding 
in  delay  to  make  an  accurate  measurement.  In  particular,  if  the  delay 


grid  size  Ax  Is  chosen  so  there  are  only  a  few  delay  bins  encompassing  the 
delay  spread  of  the  signal,  then  the  measurement  of  the  time  delay  jitter 
will  be  limited  by  the  delay  grid  size.  The  resulting  frequency  selective 
bandwidth,  which  Is  computed  from  the  time  delay  jitter,  will  be  inaccu¬ 
rate.  An  accurate  measurement  of  the  frequency  selective  bandwidth  re¬ 
quires  on  the  order  of  20  to  30  delay  bins  in  the  delay  region  which 
encompasses  95  percent  of  the  signal  power  (see  liquation  5-21).  Because 
the  delay  sample  size  Ax  Is  often  chosen  on  the  basis  of  a  channel  symbol 
period  rather  than  to  achieve  an  accurate  measurement  of  the  frequency 
selective  bandwidth,  a  comparison  of  the  signal  power  in  each  delay  bln 
with  the  ensemble  values  may  be  a  better  check  that  the  realization  has 
the  proper  delay  distribution. 

The  accuracy  with  which  the  GPSO  Is  evaluated  can  be  checked  by 
computing  the  mean  power  P_  in  the  K  -  K  -  x  grid  (Equation  5-3  or 

w  a  y 

5-35).  The* mean  power  out  of  an  antenna  should  be  equal  to  1/L^  where 
Is  the  ensemble  scattering  loss  of  the  antenna.  The  difference  in  the 
values  of  PQ  and  1 /L s  Is  equal  to  the  amount  of  delayed  signal  power  that 
Is  not  Included  in  the  delay  grid  If  the  integrals  involved  In  the  evalua¬ 
tion  of  the  mean  signal  power  In  each  grid  cell  are  done  with  sufficient 
accuracy.  If  the  mlnlmin  delay  criterion  given  by  Equation  5-21  Is  met, 
then  these  two  values  should  agree  to  within  a  few  percent. 

The  following  subsections  describe  the  algorithms  used  to 
measure  the  realization  parameters. 

5.3.1  Scattering  Loss. 

a 

The  mean  power  of  the  realizations  will  be,  in  general,  less 
than  unity  because  of  the  signal  attenuation  caused  by  the  antenna  beam 
profile.  The  difference  between  the  mean  power  of  the  realization  and 
unity  Is  equal  to  the  scattering  loss.  The  steps  for  computing  the 
scattering  loss  of  a  realization  are  as  follows: 
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1.  Compute  the  impulse  response  function  integrated  over  all 
signal  delays.  In  continuous  notation,  this  is 

hA,m(*'  =  .1  hA,m<*-''>dT  •  l5‘37> 

This  integral  is  evaluated  in  the  channel  simulator  as 

hA,m(klx)  a  11  |  hA,m<kl’<-J4,>  (5‘38> 

where  the  sum  Is  over  all  delay  bins.  The  quantity 
hA  m(kAx)  represents  the  voltage  that  would  be  received 
by  the  antenna  under  flat  fading  conditions  (e.g.  when 
the  transmitted  modulation  is  turned  off). 

2.  The  power  Pm  of  the  Impulse  response  function  is  averaged 
over  the  realization: 

N 

Pm  ■  <1/Nx>  J  |hA,m(kix)|2  •  (5-39) 

k«l 

3.  The  scattering  loss  In  dB  is  computed  for  the  mth  antenna: 

LS,m  ■  -10  log10(Pm)  .  (5-40) 

5.3.2  Frequency  Selective  Bandwidth  and  Delay  Distribution. 

The  frequency  selective  bandwidth  is  computed  by  calculating  the 
time  delay  jitter  of  the  signal.  Sufficient  delay  resolution  is  required 
to  do  this  computation  accurately.  The  realization  delay  distribution  Is 
given  by  the  mean  power  of  the  j^  delay  bln.  These  values  may  be 
compared  with  the  approximate  ensemble  values  given  by  Equation  5-3  or 
5-36.  The  algorithms  for  computing  the  frequency  selective  bandwidth  and 
the  delay  distribution  are  as  follows: 
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1.  Calculate  the  first  two  moments  of  delay: 


Pm<T">  =  d'"x>  I  -  (j^T)n|AThA  (kAx.jAx) |2  .  (5-41) 

k  j  * 

This  is  evaluated  for  n=l  and  for  n*2. 

2.  The  time  delay  jitter  is 

UT  *  [<t2>  -  <t>2]1/2  .  (5-42) 

3.  The  frequency  selective  bandwidth  of  the  realization  is 
computed  from  aT  as 

Y»  *  WJ  ■  <5-«> 

4.  Calculate  the  mean  power  In  the  jth  delay  bin: 

Pj ,m  *  (1/Mx>  1  .  (5-44) 

k 

5-3.3  X  Direction  Decorrelation  Distance  or  Time. 

The  x  direction  decorrelation  distance  or  the  decorrelation  time 
Is  computed  by  finding  the  1/e  point  of  the  autocorrelation  function  of 
the  impulse  response  function.  The  algorithm  used  to  compute  the  decor¬ 
relation  distance  or  time  Is: 

1.  Calculate  the  power  spectrum  of  the  realization.  In 
continuous  notation  this  is 

Lx^2 

W  "  |(1/Lx>  i  hA,m<x>  e*P  <-VH*  (5'45) 

-V2 

where  x  is  replaced  by  t  to  compute  the  decorrelation  time. 
The  discrete  equivalent  of  this  equation  is 
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where  ks,  the  relative  index  of  the  two  antenna  centers,  is  given  by 

k  =  (x  -xJ/ax  .  (5-52) 

s  n  m 

Here  xn  and  xffl  are  the  relative  x  coordinates  for  the  antenna  centers. 

5.*  EXAMPLES. 

In  this  section,  examples  are  presented  of  the  channel  simula¬ 
tion  techniques  outlined  In  Section  5.1.  First  a  comparison  Is  made 
between  the  ensemble  signal  parameters  and  the  values  obtained  from  random 
realizations  of  the  Impulse  response  function  at  the  output  of  an  antenna. 
These  results  show  that  the  random  realizations  do  indeed  have  signal 
parameter  values  that  are  close  to  the  ensemble  values.  Then  some  exam¬ 
ples  of  the  voltage  amplitude  at  the  output  of  a  matched  filter  to  a 
transmitted  square  pulse  are  shown  to  illustrate  some  of  the  effects  and 
techniques  that  have  been  discussed  in  this  report.  It  is  assumed  in  this 
section  that  mean  ionization  or  dispersive  effects  are  negligible  in  com¬ 
parison  to  stochastic  diffraction  effects.  It  is  further  assumed  that  the 
penetration  angle  is  zero  so  the  scattering  is  Isotropic  about  the  line- 
of-sight. 


5.4.1  Signal  Parameters  of  Random  Realizations. 

The  signal  parameters  at  the  outputs  of  the  square,  Gaussian 
antennas  without  pointing  errors  and  without  sidelobes  are  considered 
here.  The  random  realizations  of  the  impulse  response  function  were 
generated  at  the  outputs  of  two  antennas  whose  centers  are  displaced  by 
the  distance  D  where  D  Is  the  length  of  a  side  of  one  antenna.  The  rota 
tion  angle  was  chosen  to  be  45*.  The  configuration  is  similar  to  that 
depicted  in  Figure  5-1  except  that  there  are  only  two  arrays  rather  than 
the  three  shown  in  the  figure.  The  algorithms  presented  in  Section  5.3 
were  then  used  to  compute  the  realization  signal  parameters. 
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Under  these  conditions,  Equations  3-55  give  the  ensemble  scat¬ 
tering  loss,  and  f ^/f o*  The  cr0SS  correlation  of  the  impulse 

response  functions  out  of  the  two  antennas  is,  from  Equation  3-40, 

pa  =  exp  [-(D/i0)2U0/1a)2]  .  (5-53) 

Figure  5-3a  through  5-3d  present,  respectively,  the  ensemble 

(solid  lines)  and  realization  (dots)  values  of  the  scattering  loss,  fre¬ 

quency  selective  bandwidth,  x-direction  decorrelation  distance,  and  the 
voltage  cross  correlation  pa  of  the  Impulse  response  functions  at  the 
outputs  of  the  two  antennas.  Except  for  the  cross  correlation,  the  reali¬ 
zation  values  plotted  are  the  average  of  the  values  for  the  two  antennas. 

Because  of  the  finite  length  of  the  realizations,  the  agreement 

between  the  realization  values  and  the  ensemble  values  is  not  exact. 

However,  over  the  range  of  0/tQ  shown  from  0.1  to  5,  the  agreement  of  the 

scattering  loss  and  the  ratio  fA/fQ  with  the  ensemble  values  is  quite 

good.  Somewhat  more  scatter  about  the  ensemble  values  is  seen  in  the 

values  of  the  ratio  l^/lQ  and  the  cross  correlation  function.  The  ratio 

i./i  appears  to  be  biased  to  the  high  side  which  is  attributed  to  the 

a  o 

linear  Interpolation  algorithm  used  to  solve  Equation  5-50  for  the  1/e 
point  of  the  realization  autocorrelation  function. 

5.4.2  Received  Voltage. 

This  final  section  presents  examples  of  the  received  voltage  out 
of  an  integrate  and  dump  filter  that  is  matched  to  a  transmitted  square 
pulse.  These  examples  are  intended  to  illustrate  the  effects  of  frequency 
selectivity  and  antenna  filtering  on  a  transionospheric  communications 
link  and  to  illustrate  the  differences  in  the  structure  of  the  received 
signal  depending  on  whether  the  frozen-in  or  turbulent  approximation  is 
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used  to  generate  the  impulse  response  function  realizations.  The  follow¬ 
ing  calculation  will  also  illustrate  how  the  received  voltage  can  be  con¬ 
structed  from  the  impulse  response  function  realizations  in  a  digital  sim¬ 
ulation  of  the  link.  Additional  examples  for  specific  system  applications 
may  be  found  in  Bogusch,  et  al.  (1981)  and  in  Bogusch,  Guigliano,  and 
Knepp  (1983). 

The  received  time-varying,  complex  voltage  from  a  transmitted 
impulse  under  conditions  where  the  mean  ionization  effects  are  negligible 
is  given  by  the  inverse  of  Equation  2-103: 

ao 

U(w+u)  t)  =  /  h(T,t)  exp  (-iurr)dT  (5-54) 

.00 

where  h(t,t)  is  the  impulse  response  function  and  u>Q  is  the  carrier  fre- 
auency.  The  quantity  U( CJ+a)0 » t )  then  represents  the  voltage  at  a  relative 
frequency  w  out  of  an  antenna  at  time  t.  For  a  transmitted  square  pulse 
with  a  chip  duration  T  ,  the  voltage  out  of  the  matched  filter  at  time  t 

v 

can  be  written  as 

00 

e(T,t)  =  (1/2*)  /  M(u))U(u)+(i>  ,t)  exp  (iu>T)dw  (5-55) 

«00 

where  t  is  the  relative  time  delay  of  the  matched  filter  and  where 

M(w)  =  T  sin2(»T  /2)/(mT  /Z)2  (5-56) 

is  the  combined  spectrum  of  the  square  pulse  and  the  matched  filter. 

In  an  actual  simulation,  these  two  integrals  are  performed  at 
baseband  (^=0)  and  must  be  performed  discretely.  If  the  impulse  response 
function  is  generated  with  Nq  delay  samples  of  size  At,  then  the  dis¬ 
crete  frequency  domain  voltage  at  time  t  is 
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u 

U(kAu),t)  =  l  Aih( j  At  ,  t )  9X0  [  -i  k Aw( t^+jAt  )  ]  (5-57) 

j=0 

where  Aw  =  2*/Ndat.  Here  the  delay  ts  defines  the  start  of  the  delay  grid 
upon  which  the  impulse  response  function  is  generated  (Equation  5-1). 
However  this  delay  may  be  ignored  at  this  point  because  it  will  be 
accounted  for  when  the  received  voltage  is  constructed.  The  range  of  the 
index  k  in  this  equation  is  from  0  to  Nq-1  so  the  range  of  frequencies 
represented  is  from  0  to  (Ng-l)Au.  The  zero  frequency  of  U(kAw,t)  must 
correspond  to  the  smallest  frequency  of  M(u)  within  the  receiver  band¬ 
width.  The  discrete  received  voltage  is  then  given  by 

N°“i  sin2 [ ( k-Nn/2) AwT  /2] 

e( x,t)  =  ( AwT  /2it)  l  - -  U(kAw,t)  exp  (ikAw(T-T')]  . 

C  k=0  [ ( k-NQ/2) AwTc/2] 2  (5-58) 

If  At  is  chosen  to  be  T./2,  then  AwTc/2  =  2n/ND  and  e(x,t)  represents  a 
signal  that  is  band-limited  to  the  frequency  range  -l/T.  to  -t-i/T..  Note 
also  that  e(t,t)  is  unambiguous  in  delay  only  over  the  interval  from  0  to 
(Nq-I)At.  The  arbitrary  starting  delay  t1  is  included  in  Eauation  5-49  to 
account  for  the  delay  t$  which  may  have  been  ignored  in  Equation  5-48  and 
to  ensure  that  e(x,t)  is  not  aliased  in  delay. 

In  the  examples  that  follow,  the  random  realizations  of  the 
impulse  response  function  were  generated  with  a  delay  sample  size  of  Tc / 2 
and  with  32  delay  samples.  The  chip  rate  Rc  =  1/TC  was  set  at  1  MHz. 
However,  the  frequency  selective  effects  depend  only  on  the  ratio  of  the 
frequency  selective  bandwidth  to  the  chip  rate  fQ/R  .  Then  the  voltage 
out  of  a  matched  filter  was  calculated  from  the  impulse  response  functions 
usinq  Equations  5-48  and  5-49.  The  amplitude  of  this  voltage  is  plotted 
as  a  function  of  the  delay  of  the  matched  filter  and  as  a  function  of 
time . 
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The  effect  of  frequency  selectivity  on  the  received  voltage 
amplitude  is  shown  in  the  next  set  of  figures.  The  impulse  response 
functions  used  for  these  signals  were  generated  using  the  frozen- in 
approximation  and  an  omnidirectional  antenna.  Delay  in  chips  is  plotted 
on  the  abscissa  of  the  figures  and  time  or  antenna  position  is  directed 
into  the  figure.  The  total  time  duration  of  the  plots  Is  10  t0. 

The  signal  in  Figure  5-4a  where  f  /R  *  1.0  clearly  shows  the 
effect  of  signal  fading  as  the  peak  amplitude  rises  and  falls  with  time. 
The  matched  filter  output  in  the  figure  is  somewhat  rounded  rather  than 
being  a  triangle  because  the  signal  is  band-limited  to  the  frequency  range 
between  the  first  nulls  of  the  transmitted  spectrum.  Some  minor  distor¬ 
tion  of  the  output  waveform  is  seen  but  for  the  most  part  the  signal  is 
contained  within  the  period  of  one  transmitted  symbol.  For  this  channel 
the  fading  is  nearly  flat  which  means  that  all  frequency  components  within 
*  the  signal  bandwidth  propagate  essentially  the  same  way  throuqh  the  dis¬ 
turbed  ionosphere  and  there  is  very  little  time  delay  spreading  of  the 
received  signal . 

Figure  5-4b  shows  a  signal  generated  with  fQ/R  =  0.5.  The 
impulse  response  function  of  this  signal  was  generated  using  the  same  set 
of  random  numbers  and  the  same  angular-delay  grid  that  was  used  for  the 
Impulse  response  function  of  the  signal  in  the  previous  figure.  The  only 
difference  in  the  two  impulse  response  functions  is  the  distribution  of 
signal  energy  within  the  delay  bins.  Thus  there  is  strong  correlation  in 
the  signal  structures  seen  in  the  two  figures.  For  this  smaller  value  of 
the  frequency  selective  bandwidth,  more  of  the  signal  energy  is  arriving 
at  longer  delays.  There  is  one  time  about  three  quarters  of  the  way  into 
the  figure  where  there  are  two  distinct  peaks  in  the  matched  filter  output 
amplitude.  It  is  these  sort  of  structures  that  can  cause  delay  tracking 
algorithms  in  receivers  to  lose  lock  and  that  cause  intersymbol  inter¬ 
ference  which  can  degrade  demodulation  performance. 
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The  matched  filter  output  amplitudes  for  f  /R£  equal  to  0.2  and 
0.1  are  shown  in  Figures  5-4c  and  5-4d  respectively.  The  signals  in  these 
figures  are  much  more  distorted  than  the  signals  in  the  previous  two 
figures  with  the  signal  energy  being  spread  out  over  multiple  chips.  This 
amount  of  multipath  delay  spread  produces  severe  intersymbol  interference 
and  reduces  the  amplitude  out  of  the  matched  filter  at  a  given  delay. 

Another  effect  that  Is  evident  In  these  signals  is  that  the 
signal  received  at  long  delays  varies  more  rapidly  in  time  than  the  signal 
received  at  shorter  delays.  This  is  a  result  of  the  frozen- in  approxima¬ 
tion.  Figure  5-5  shows  a  signal  with  f  /R  *  0.1  whose  impulse  response 
function  was  generated  using  the  turbulent  approximation.  Again  10 
decorrelation  times  of  the  signal  are  plotted.  It  should  be  remembered 
when  comparing  the  signals  in  Figure  5-4d  and  5-5  that  the  Doppler 
spectrum  has  a  Gaussian  form  for  the  frozen-in  approximation  and  an  f‘4 
form  for  the  turbulent  approximation.  It  is  this  fact  that  accounts  for 
the  more  splkey  appearance  of  the  turbulent  approximation  signal.  Also 
the  impulse  response  function  generation  techniques  are  sufficiently 
different  so  there  is  no  correlation  between  features  seen  in  the  two 
figures.  The  thing  to  note  about  the  turbulent  approximation  signal  is 
that  the  features  at  long  and  short  delays  have  the  same  decorrelation 
time. 

Finally,  the  effect  of  an  antenna  on  the  received  signal  ampli¬ 
tude  is  shown  in  the  last  set  of  figures.  The  Impulse  response  functions 
used  for  these  signals  were  generated  using  the  frozen-in  approximation 
and  a  square  Gaussian  antenna.  As  is  the  case  for  all  signals  plotted  in 
this  section,  the  scattering  is  isotropic  about  the  1 ine-of-sight.  All  of 
these  signals  have  the  ratio  f  /R„  set  at  0.1  while  the  ratio  D/JL  is 
varied.  Again  the  signal  amplitude  out  of  the  matched  filter  is  plotted 
versus  delay  with  time  directed  into  the  figure.  The  total  time  duration 
of  the  plots  is  10  -c0> 
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DELAY  (CHIPS) 

Figure  5-5.  Matched  filter  output  amplitude  for  a  signal  generated  using 
the  turbulent  approximation  ( f 0 / Rc  =  0.1). 

Figure  5-6a  shows  the  matched  filter  output  amplitude  for  a 
Gaussian  square  antenna  wl  ;.h  DA0  =  0.5.  For  this  antenna  and  isotropic 
scattering,  the  ensemble  scattering  loss  Is  0.4  dB.  This  signal  Is  not 
Identical  to  that  shown  in  Figure  5-4d,  although  similarities  can  be  seen, 
because  of  the  antenna  filtering  and  because  the  angular-delay  grid  sizes 
depend  on  the  ensemble  signal  parameters  at  the  antenna  output.  Figures 
5-6b  through  5-6d  show  respectively  the  matched  filter  output  amplitude 
for  D/*o  equal  to  1.0  (ensemble  scattering  loss  equal  to  1.3  dB),  2.0 
(ensemble  scattering  loss  equal  to  3.9  dB),  and  5.0  (ensemble  scattering 
loss  equal  to  10.0  dB) . 

The  effects  of  an  antenna  on  the  signal  parameters  are  readily 
apparent  from  these  figures  as  the  ratio  D/i0  is  reduced.  The  signal 
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energy  arriving  at  longer  delays  Is  preferential^  suppressed  by  the 
antenna.  For  D/A  equal  to  5.0,  the  output  signal  looks  almost  flat  with 
very  little  distortion  of  the  matched  filter  output.  However,  the  output 
signal  anplltude  Is  10  dB  smaller  on  the  average  than  it  was  without 
antenna  filtering.  Also,  the  decorrelation  time  or  distance  of  the 
signals  is  seen  to  Increase  as  the  antenna  filtering  Increases. 
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APPENDIX  A 

PHASE  VARIANCE  OUE  TO  ELECTRON  DENSITY  FLUCTUATIONS 


A  relationship  between  the  phase  variance  imparted  on  the  wave 
as  it  propagates  through  the  ionization  layer  and  the  electron  density 
fluctuations  will  be  derived  in  this  appendix.  This  relationship  is  given 
by  Equation  2-51  which  was  derived  using  the  Markov  approximation. 

However,  It  will  be  shown  here  that  the  relationship  requires  only  that 
the  layer  thickness  be  large  compared  to  the  decorrelation  distance  of  the 
electron  density  fluctuations  along  the  1 ine-of-sight. 

It  was  shown  in  Section  2  that  the  total  phase  change  of  the 
wave  as  it  propagates  through  the  ionization  layer  is 

=  rex<ne>  /  5(p,z,t)dz  .  (A-l) 

The  autocorrelation  function  of  the  phase  fluctuations  is  then 

<4»(p,tH(p't’)>  =  (reX<ne>)2  J  dz  /  dz'  <5(p,z,tH(p\z\t')>  (A-2) 

O  0 

where  t  is  the  thickness  of  the  scattering  layer.  For  spatially  and 
temporally  stationary  random  electron  density  fluctuations,  the  expecta¬ 
tion  must  be  a  function  of  the  differences  Ap,  z-z'  and  At  only.  Denoting 
the  autocorrelation  of  5(p,z,t)  by  B^(Ap,z-z' ,At)  and  the  autocorrelation 
of  4> ( p , t )  by  B^( Ap, At) ,  Equation  A-2  becomes 

♦  ,  L  L 

B  (Ap.At)  =  (r  Kn  >)  /  dz  /  dz’  B£(Ap,z-z ' , At)  .  (A-3) 

?  oo 
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This  double  integral  may  be  reduced  to  a  single  integral  by  changing  the 
order  of  integration  with  the  result 


If  the  correlation  distance  of  B^Ap.z.At)  along  the  z  direction 
is  small  compared  to  L,  then  B^(Ap,z,At)  will  become  small  before  1 2 j /L 
approaches  unity  in  the  integral  and  the  1 2 1 /L  term  may  be  ignored.  The 
limits  of  the  integral  may  then  be  set  to  ±  ®  and  the  integral  reduces  to 

m 

B^(Ap ,  At)  =  (rex<ne>)2  L  /  B£(&p,z, At)dz  .  ( A-5 ) 

The  remaining  integral  is  denoted  by  A(Ap,At)  so  the  autocorrelation  of 
the  phase  fluctuations  is 

B . ( Ap, At)  =  (r  X<n  >)2  L  A(  Ap, At)  (A-6) 

C  C 

which  is  the  same  as  the  final  expression  in  Equation  2-51. 
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APPENDIX  8 

SIGNAL  PARAMETERS  FOR  ELECTRON  DENSITY  FLUCTUATIONS 


In  this  appendix,  the  expansion  coefficients  A0  and  A2  are 
calculated  using  the  quadratic  approximation  of  the  correlation  function 
A(P(j)  of  the  electron  density  fluctuations  and  usina  the  delta  layer 
approximation.  From  these  coefficients,  the  phase  variance,  decorrelation 
distance,  and  the  coherence  bandwidth  of  the  signal  incident  on  the  plane 
of  the  receiver  are  written  in  terms  of  physical  parameters.  However, 
these  signal  parameters  are  computed  from  a  disturbed  ionosphere  model 
using  the  more  general  formalism  of  Wittwer  (1979,  1980)  which  accounts 
for  the  finite  thickness  of  the  scattering  region  and  other  complicating 
effects.  The  purpose  of  this  appendix  is  only  to  illustrate  the  depend¬ 
ence  of  the  signal  parameters  on  geometrical  and  electron  density  fluctua¬ 
tion  parameters. 


A  power- law  form  of  the  PSD  for  the  three  dimensional  electron 
density  fluctuations  is  assumed: 


where 


8*3/2lrlsLt<An!>r(n) 

<ne>2r(n-3/2)(l+K-L-K)n 


( B-l ) 


rr 

j  ^ 
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The  scales  Lr,  l$,  and  L  ,  Lz,  L^z  are  defined  in  Section  2.6  in 
terms  of  the  outer  scale  l0,  the  penetration  angle  »!>,  and  the  axial 
ratio  q.  The  root  mean  square  value  of  the  electron  density  fluctuations 
is  <An2>.  For  K*l*K  >  1,  is  proportional  to  «'3n,  Thus  a  K*4  PSD  for 
the  three  dimensional  electron  density  fluctuations  corresponds  to  the 
n  =  2  case. 


The  correlation  function  A(pd)  under  the  delta  layer  approxi¬ 
mation  is  given  by  Equation  2-70  which  is  reproduced  here: 

00 

A(Pd)  *  (2u)-2  //  exp(ii^«Pd)$c(i<j^Kz=0)d2i<j^  .  (2-70) 


Using  Equations  B-l  and  R-3  and  performing  the  angular  integral,  Equation 
B-3  reduces  to 

♦  4/it  q6L  T( n)<An2>  *  Jo^Pd*L  I  *V  / 

A(Pd)  «  - - -  / -  u  du 

T(n-3/2)<ne>2  o  (l+u2)n 


■  ^.;<LQ(n:1-<.Anl>  (Pd-lVPd)(n'1)/2  Kn_1[(Pd-L,-Pd)1/2]  (B-3) 

2n‘3r(n-3/2)<ne>2  -L 


where  JQ  is  the  Bessel  function  of  order  0,  is  the  modified  Bessel 
function  of  order  n-1,  1^  is  a  2X2  matrix  containing  the  x-y  components 
of  L,  and  6  is  defined  in  Equation  2-69. 


For  all  values  of  n  except  n  =  2,  A(pd)  can  be  expanded  in  a 
power  series  of  the  form 

HoJ  •  A0((1-A2[{x/Lx)2  +  (y/uy)2]m/2)  (B-4) 
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where  m  =  min(2,2n-2).  For  n  =  2,  A2  does  not  exist  unless  an  inner  scale 
li  is  imposed.  This  is  accomplished  by  truncating  the  integral  over  u 
in  Equation  B-3  at  a  cutoff  u  =  where  e  is  chosen  to  be  1^/L  .  For 
values  of  n  greater  than  2,  the  J0  Bessell  function  in  the  integrand  of 
Equation  B-3  can  be  expanded  and  the  resulting  series  can  be  integrated 
term-by-term.  The  first  two  terms  of  the  expansion  give 


2A  q6L  r(n)<Anf> 


A  = 
o 


r(n-3/2)<ne>' 


[l-l/U+e'2)0"1] 


(3-5) 


and 


A2  = 


.  ( l+e”2 )n~1  -  [l+(n-l)e'2] 


(B-6) 


4(n-2) ( ( 1+e"2) n-1  -1. ] 

In  the  limit  that  n  =  2,  the  A2  coefficient  becomes 


A?  = 


( 1+e"2 ) 1 n ( l+e~2 )  -  e'2 


4e 


-2 


(B-7) 


This  expression  can  be  further  reduced  in  the  limit  that  li  «L0  to 

A2  -  ln(L0/l.)/2  .  (B-8) 

Now  the  phase  variance,  the  decorrelation  distance,  and  the 
coherence  bandwidth  can  be  written  in  terms  of  geometrical  parameters  and 
electron  density  fluctuation  parameters.  Using  Equation  2-52,  the  phase 
variance  due  to  the  structured  ionization  is 


-  2q«L0b6(reA )2<An*>  (B-9) 

where  is  the  thickness  of  the  delta  layer.  The  decorrelation  distance 
and  the  coherence  bandwidth  are  defined  in  Equations  2-78  and  2-79 
respectively.  Using  these  definitions. 
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•*  * 


(B-10) 


and 


«  . 
coh 


;’,V<Vzr)l-o 

lztzrln(Lo/’tK 


(B-ll) 


These  equations  are  only  valid  for  the  delta  layer  approxima¬ 
tion,  for  the  quadratic  phase  structure  approximation  and  for  a  K“4  three 
dimensional  electron  density  fluctuation  PSD.  It  can  be  seen  from  the 
equations  that  the  values  of  i  and  <ocoh  are  only  weakly  dependent  on  the 
inner  scale  1  j . 
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